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Abstract 

We study the hadron spectra in nearly central A+A collisions at RHIC and LHC in a 
broad transverse momentum range. We cover the low-pr spectra using longitudinally 
boost-invariant hydrodynamics with initial energy and net-baryon number densities 
from the perturbative QCD (pQCD)+saturation model. Build-up of the transverse 
flow and sensitivity of the spectra to a single decoupling temperature T dec are studied. 
Comparison with RHIC data at v /J/v r /v = 130 and 200 GeV suggests a rather high 
value T dec = 150 MeV. The high-p T spectra are computed using factorized pQCD cross 
sections, nuclear parton distributions, fragmentation functions, and describing partonic 
energy loss in the quark-gluon plasma by quenching weights. Overall normalization is 
fixed on the basis of p+p(p) data and the strength of energy loss is determined from 
RHIC Au+Au data. Uncertainties are discussed. With constraints from RHIC data, 
we predict the pt spectra of hadrons in 5 % most central Pb+Pb collisions at the LHC 
energy ^snn = 5500 GeV. Due to the closed framework for primary production, we can 
also predict the net-baryon number at midrapidity, as well as the strength of partonic 
energy losses at the LHC. Both at the LHC and RHIC, we recognize a rather narrow 
crossover region in the pr spectra, where the hydrodynamic and pQCD fragmentation 
components become of equal size. We argue that in this crossover region the two 
contributions are to a good approximation mutually independent. In particular, our 
results suggest a wider p^-region of applicability for hydrodynamical models at the 
LHC than at RHIC. 
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1 Introduction 



Transverse momentum spectra of hadrons in ultrarelativistic heavy ion collisions offer 
invaluable input for the studies of QCD matter properties and QCD dynamics. At high 
collision energies, currently at BNL-RHIC and in the near future at the CERN-LHC, 
the parton and hadron production rates are significant over a wide transverse momen- 
tum range. The high- and low-pT parts of the transverse momentum spectra provide 
complementary information on the QCD matter. As reviewed below, the hadron spec- 
tra measured in Au+Au collisions at RHIC show systematic evidence on formation of 
elementary-particle matter and provide motivation for hydrodynamical description of 
A+A collisions. 

Large energy densities. First, average initial energy densities can be estimated 
on the basis of measured transverse energy |T] originating from the low-py region of 
the hadron spectra. Even the conservative Bjorken estimate [2|, e = ^^/(ttR^t), 
computed using dEr/dr] = 606 . . . 620 GeV as measured at r\ ~ in ... 5 % central 
Au+Au collisions at y/s^N = 200 GeV at RHIC, implies an average local energy 
density e ~ 5 GeV/fm 3 at r ~ 1 fm, a typical hadronization time in a p+p collision 
in vacuum. Already this is well above the energy densities where hadrons could form. 
Since the formation time of such a system at RHIC is most likely considerably shorter 
than 1 fm, the average initial energy densities can be expected to be clearly, say, even 
a factor five, higher. Furthermore, if there is pressure in the system, the estimates 
of energy density become even larger, since dEx/dy is not conserved due to the pdV 
work by pressure in longitudinal expansion. According to the hydrodynamic studies, 
the initial-to-final transverse energy degradation can be as much as a factor three [HJ, 
increasing the initial energy density estimates accordingly. 

Pressure. Second, clear azimuthal asymmetry in the transverse spectra has been 
observed [3] in non-central collisions. This phenomenon is possible only if sufficiently 
frequent secondary collisions take place to build up pressure. Furthermore, hydrody- 
namic studies imply that the pressure must form early enough (in less than 1 fm) 
for the spatial azimuthal asymmetry in the primary production to be transferred into 
momentum distributions 5\. Complete thermalization, however, is not necessary for 
generating pressure in the system. Both observations, the large energy densities and 
the early existence of pressure, considered together suggest that a collectively behaving 
partonic matter of gluons, quarks and antiquarks, is produced. 

Partonic energy loss. Third, the high-pr hadron spectra in Au+Au collisions at 
Pt ft 4 . . . 5 GeV have been observed to be suppressed by a factor ~ 4 ... 5 relative to 
P+P S3 El El EE EDI EI]- This energy loss, jet quenching, is one of the predicted QGP 
signals [T2] . The fact that no such suppression is discovered in d+Au collisions [T5] . 
excludes the possibility that the primary QCD collision mechanism for producing high- 
er partons in nuclear collisions would be very different from that in p+p. Also the 
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back-to-back correlation of the high-p^ hadrons, present in the p+p collisions, has been 
observed to disappear in central Au+Au [H]. These results, considered together with 
the first two, very strongly suggest that the origin of the observed high-p T suppressions 
in Au+Au is the produced partonic medium: before hadronization the hard partons 
lose energy in penetrating through the dense parton matter. 

Hydro for low-pT and pQCD for high-p^ spectra. The interpretation of the 
HBT radii extracted from the particle correlation measurements ^3] is still a puzzle to 
be solved for a better understanding of the produced collective system. In spite of this 
caveat, the RHIC results indicate that a dense strongly interacting partonic system, 
where mean free paths are much smaller than the size of the system, is produced in 
nonperipheral A+A collisions. It is then reasonable to expect that hadron spectra 
at pr<2...3 GeV, the region which dominates the total multiplicity and transverse 
energy, can be obtained by treating the expansion of produced matter quantitatively 
in terms of relativistic hydrodynamics JH]- Several pieces of input are necessary for 
such calculation: initial state of the evolution must be given in detail, an Equation of 
State (EoS) is needed to close the system of hydrodynamic equations and finally the 
collective description of flowing matter must be turned into particle spectra. 

In the hydrodynamic approach also the originally produced high-pr partons are as- 
sumed to be absorbed into the thermalized system and all final state hadrons are 
thermally emitted from the decoupling surface. Especially, no remnants of pQCD 
power-law tails are left on top of the steeply falling thermal spectra of hadrons. In 
the spectra of hadrons measured in Au+Au collisions at RHIC, the powerlaw like 
tails do appear at large transverse momenta, pr>5..6 GeV, but since the high-pr 
partons carry only a small fraction of the total (transverse) energy, the assumption 
that all partons thermalize does not affect the treatment of the thermal part apprecia- 
bly. For the high-p^ region, a more relevant baseline calculation is hadron production 
without rescatterings, where all the produced high-p^ partons cascade in vacuum into 
non-interacting partonic showers which then independently hadronize. Such reference 
spectra can be computed for A+A collisions using perturbative QCD (pQCD) in the 
same way as in p+p and p+p collisions. Energy losses can also be incorporated into 
these baseline spectra [T2] . 

In this paper we aim at a comprehensive description of hadron spectra in central 
A+A collisions, comparing the results with RHIC data and giving detailed predictions 
for the 5 % most central Pb+Pb collisions at the LHC. On one hand, we compute the 
spectra within the framework of hydrodynamics, applying initial conditions predicted 
by pQCD and saturation [T71 E], and also including decays of produced resonances 
after decoupling. On the other hand, we compute the high-p^ reference spectra in the 
framework of collinearly factorized pQCD cross sections, nuclear parton distributions 
(nPDFs) and fragmentation functions JH] • Following Ref . JH] , energy losses of high- 
Pt partons are incorporated in terms of quenching weights [201 • A comparison of the 
computed low- and high-pr spectra with RHIC data is presented and, in particular, 
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these spectra are predicted for the 5% most central Pb+Pb collisions at the LHC. 
We also determine and discuss the theoretically challenging crossover region where the 
hydrodynamic and the pQCD spectra are comparable. 

The rest of this paper is arranged as follows: In Sec. 01 we first briefly review the 
saturation models from which the initial state for hydrodynamics can be computed at 
collider energies. After this, we discuss in more detail the pQCD+saturation initial 
conditions. In Sec. 01 we define the hydrodynamic framework. In particular, we discuss 
the fast thermalization assumption, the EoS, and collective flow phenomena. Section 0] 
specifies the pQCD computation of the baseline high-p^ spectra without and with 
energy losses. The results are presented in Sec. 01 and conclusions in Sec. 01 

2 Initial state for hydrodynamics 

A key input for hydrodynamical studies is the initial state of evolution. At the SPS 
energies, hadron multiplicities and spectra have been used to constrain the initial state, 
the EoS, flow and decoupling [2TJ02]- However, it turned out that different combina- 
tions of initial states and EoS can lead to similar results for the final spectra. Hence, 
the spectra could obviously be more efficiently used to extract information on the prop- 
erties of QCD matter and its spacetime evolution if the initial state could be extracted 
from elsewhere. This can be expected to be the case at collider energies as perturba- 
tively computable parton production may dominate initial parton production in A+A 
collisions at sufficiently high collider energies j2Hl EH 123 EHJ Due to the diverging 
partonic cross sections pQCD alone, however, does not predict the amount of trans- 
verse energy deposited into the central rapidity region. A further element, like gluon 
saturation, is needed. We briefly comment on the saturation models which provide a 
basis for estimating the initial conditions. 

2.1 A brief review of saturation models 

2.1.1 The general idea in saturation - evidence in DIS 

Gluon saturation was introduced more than two decades ago for p+p collisions [27| as a 
property of the wave functions of the incoming hadrons which makes the primary gluon 
production finite and preserves unitarity. At sufficiently small momentum fractions x 
and small scales Q 2 , the evolution dynamics of gluon distributions (BFKL [2E] or 
DGLAP [21]) is expected to become dominated by gluon fusions rather than splittings, 
inhibiting the growth of the gluon distributions towards smaller x or higher Q 2 and 
leading to saturation. The idea has been successfully tested in the case of the free 
proton: geometric scaling of the structure function F$ at small x and small Q 2 predicted 
by the dipole picture [30] has been observed. Also the leading-order (LO) DGLAP 
fits to the HERA data in the region x ~ 10~ 5 and Q 2 ~ 1.5 GeV 2 improve after 
including nonlinearities [23 01] induced by the gluon fusions [321 ■ Extrapolation to 
nuclei, however, is more model-dependent, see e.g. [331 Oil- 
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2.1.2 Saturation in A-\-A - semianalytical models 

In the context of ultrarelativistic heavy ion collisions, the idea of gluon saturation in 
the wave functions of the colliding nuclei governing and regulating the final state par- 
ticle production was discussed first in [2H]. It was suggested that particle production 
in high energy A+A collisions could be computed in perturbative QCD (pQCD) as 
the saturation scale Q sat oc A 1 / 6 becomes large, Q sat 3> Aq CD . More recently, initial 
gluon production in A+A collisions was computed semi-analytically in Refs. [33] HE] 
by applying kx factorization to 2 — > 1 processes. In this approach, the shape of final 
state rapidity distributions is taken to be directly that of initially produced gluons, 
and the overall normalization is fixed on the basis of RHIC data at one energy. Good 
fits to the multiplicities and pseudorapidity distributions at RHIC have been obtained 
[TTH| and predictions for the LHC have been made 36J. It is, however, useful to recall 
the theoretical uncertainties of this approach. Widening of the rapidity distributions 
due to secondary collisions, as demonstrated by longitudinally boost-noninvariant hy- 
drodynamics [37J EE1 122] j is neglected. This phenomenon affects off-central rapidities 
more than the central ones. There are also uncertainties in the determination of the 
saturation region for nuclei (see also as well as in the specific form of the uninte- 
grated nuclear gluon distributions in the saturation region and also at large x (see the 
discussion in |3T?j). 

2.1.3 Colour glass condensate model for A-\-A - a lattice approach 

The idea of saturation in ultrarelativistic heavy ion collisions was exploited also in the 
colour-glass condensate (CGC) model [10], where the small-x gluons are described by 
classical gauge fields. Over the recent years this model has experienced a very active 
stage of development, for recent reviews see e.g. [33]. With boost-invariance the CGC 
model becomes an effective 3-dimensional theory. Primary gluon production can then 
be computed nonperturbatively by assuming a Gaussian distribution of random colour 
charges and solving the classical Yang- Mills equations on a lattice. For the SU(3) 
case, this was first done in [42J. The too large transverse momentum relative to the 
multiplicity obtained from the calculation was, however, a puzzle (see e.g. (32]) which 
was solved only recently in 03] . The lattice approach, however, does not provide a value 
for the saturation scale, and thus no absolute prediction for the initial state, either. As 
discussed in 0133], the predicted initial state of matter can range between (i) a dilute 
system of gluons, which only fragment into more gluons but do not interact with each 
other (as in the scenario of Ref. (HH! discussed above), and, (ii) a very dense system 
where gluons interact vividly and where the initially produced multiplicity is essentially 
the final multiplicity but transverse energy is clearly larger than the final Et- The latter 
case is equivalent to what is predicted in the pQCD+saturation+hydrodynamics model 
[T7| |3] discussed below. The nonperturbative CGC model [33] or the semi-analytical 
Ut factorization model (35] EH] can be used to predict the initial state at the LHC if 
the saturation scales and unintegrated gluon densities can be extracted and the form 
of colour density fluctuations constrained on the basis of deep inelastic data for 
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2.1.4 pQCD+final state saturation model in A-\-A 

In the pQCD+saturation model a dense initial state with large transverse energy 
is predicted by requiring saturation among the produced gluons [T7| EH] m the same 
spirit as was originally suggested [23 for the initial state. Parton production in 
a midrapidity unit Ay in a central A+A collision is computed by applying collinearly 
factorized pQCD cross sections with nuclear parton distributions obtained from a global 
DGLAP fit jUj to nuclear deep inelastic scattering and Drell-Yan data. The parton 
production is allowed to continue from large pt partonic jets, where pQCD is known to 
work, down to a semi-hard transverse momentum cut-off p$ ^> Aqcd- If the production 
is abundant enough the cut-off po can be defined from saturation as follows: let the 
number of produced partons (gluons clearly dominate) at px > Po in a central (b = 
0) collision be Naa(0,pq, Ay, y/s) and the effective transverse area occupied by each 
parton to be proportional to ir jp\. At saturation the partons fill the whole transverse 
overlap area nR\, and the effective saturation scale po = p sa .t is a solution of a saturation 
criterion 

N AA (0, po, Ay = 1, ^i) ■ rr/p 2 = nR 2 A . (1) 

We have taken the possible constant of proportionality in front of n jp\, which might 
contain e.g. a s or group theoretical colour factors, equal to one. 

Gluons with px < p sa ,t form at later times, r > l/p sat , and they fuse together with 
gluons formed earlier. In this sense, the number of produced partons saturates. The 
number of partons at saturation, NAA(p s &t), gives a fair estimate of the total multiplicity 
of produced partons, provided that p sat 3> Aqcd- The formation time of the dominant 
part of the system is then r ~ l/p S at- With the saturation condition of Eq. (JTJ), 
the computation of the initial state is closed and can be done at any energy. Total 
transverse energy deposit and initial energy density at r can be computed including 
the next-to- leading order contributions in minijet production |2H]. Benefits over the 
other approaches discussed above, are that chemical decomposition and especially the 
net-baryon number content of the initial state are also obtained [49J. 

The pQCD + (final state) saturation approach obviously does not suggest a micro- 
scopic dynamical mechanism for saturation, although the general idea of produced glu- 
ons to act as a medium for the to-be-produced gluons, is similar to the self-screening 
suggested in [50). The pQCD+saturation model, rather, serves as a useful effective 
tool in estimating the magnitude of initial parton production in simplest situations 
where only one scale dominates the production, as in central and nearly central col- 
lisions at midrapidities. In spite of the possibly large theoretical uncertainties, the 
pQCD+saturation model, combined with hydrodynamics, i.e. isentropic spacetime 
evolution, correctly predicted the multiplicities measured in nearly central Au+Au col- 
lisions at RHIC at y/s' = 56, 130 and 200 GeV [HI El EH- 

Also the low-pr spectra of pions, kaons and protons measured in 5% central Au+Au 
collisions at y/s = 130 GeV at RHIC are simultaneously reproduced surprisingly well 
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by pQCD+saturation+hydrodynamics, provided that the hydrodynamic evo- 
lution decouples at rather high temperature, T dec rs 150 MeV. Such a high decoupling 
temperature has also been suggested in a fireball parametrization [53]. Motivated by 
the success of the pQCD+saturation+hydrodynamics model, we here extend the study 
of Ref. to y/s = 200 GeV at RHIC and ^/~s = 5500 GeV at the LHC as well. 

2.2 Initial densities from pQCD + saturation 

According to collinear factorization and LO pQCD, the inclusive cross section of two- 
parton production in A+B collisions is given by 

^x 1 fi / A(xi,Q 2 )x 2 f j /B(x2,Q 2 )- I {s,t,u), (2) 



dp 2 T d yi dy 2 ^-^v-*,- v~-^ > dl 

where and fym are the number densities of partons % and j in the beams A and B, 
and Xi 7 2 are the fractional momenta of i and j. The hatted symbols are the Mandelstam 
variables of the partonic 2^2 subprocesses, whose cross section a oc a 2 . The partons 
k and I are produced back-to-back in transverse plane, with transverse momentum pt 
for each, and rapidities y\ and y 2 . We choose the factorization and renormalization 
scales to be identical, and set Q = px- In certain cases, such as jet production [33] and 
minijet transverse energy production [3K], the computation can be rigorously extended 
to next-to-leading order (NLO) pQCD. In these K factor which simulates the 

NLO contributions, can be defined and determined. 

Consider now only partons above a minimum scale po ^> Aqcd where the pQCD 
computation is still trustworthy. Using Eq. (j2J), the number (N?) and transverse energy 
(E T ) distributions for partons of flavour / falling into a rapidity interval Ay can be 
defined as E 

da 



dXf 



da AB ^ kl + x i 



6 (x> - p'r""' [6, t e( yi e Ay) + S„9(y 2 e Ay)]) 0(p T > p„) 



where the power n(X^) = when X? = , and n(X^) = 1 when X? = E' T . The first 
momenta of these distributions become 

cr(X f ) AyiP0 = [ dp 2 T dy l dy 2 Y] x 1 f i/A (x 1 , Q 2 )x 2 f j/B (x 2 , Q 2 ) : - j (4) 

J ij 1 + Okl 

(ki) 

da ij ^ kl ~ da ij ^ hl 
Sfk-jr {s,t,u)+5fi—z (s,u,t) 
dt dt 



p n T {Xf) 9( yi eAy)9(p T >p ) 



For an A+B collision at an impact parameter b, the average number of partons 
produced with y e Ay and px > Po can then be obtained as 

N f AB (b, Po , Ay, y/s) = T AB (b)a(N f ) 

Ay,p (5) 
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and their transverse energy as 

(6) 

The nuclear overlap function, 

T AB (b) = J d 2 r T A (h - r)T B (r), (7) 

where r = (x, y), is defined as usual, through the thickness function Ta(t) = J dzpA^z, r), 
where pa(z,t) is the nuclear density. We shall use the spherically symmetric Woods- 
Saxon nuclear density [SSI- Below, the flavour-summed quantities (N and Et) appear 
without the index /. The net-baryon number carried into the rapidity acceptance Ay 
by the px > Po minijets is then one third of the produced net quark number J3H], 

N*z*(b,p ,Ay,V~s) = TAB(b)W[*{N*) AyjPo -a(N*) Ay>po } (8) 

6 q 

= T AB (h)a{N B -*)^ P0 . 

As discussed in the previous section, the dominant transverse scale in primary pro- 
duction, po = Psat, in central rapidities of ultrarelativistic A+A collisions is obtained 
by solving Eq. Since the predictions of multiplicities at RHIC agreed with the 
measurements, we do not further tune the initial state computation from that in Ref. 
[3], i.e. we use the GRV94 parton densities [57] together with nuclear effects from the 
EKS98 parametrization [UJ, and K factors 1.6 (2.3) for RHIC (LHC) based on the 
NLO analysis of minijet Et production jlHl- Note that the i^-factors have been shown 
to vary with the PDFs in a manner that the changes in (j{Et) remain small even if 
newer PDFs were applied |15] . 

Table ^ shows the saturation scales, formation times and the pQCD quantities com- 
puted at po = p sat for RHIC and LHC in central and nearly central A+A collisions. 
The centrality selection is simulated as in Ref. by considering a central collision 
of an effective nucleus A e g < A: first, e.g. for 5 % most central collisions, we find a 
maximum impact parameter by taking 5 % of the A+A total cross section in a (op- 
tical) Glauber formulation. Then, we determine A e ^ by requiring that the number of 
participants in a central A efi +A eS collision equals the average number of participants 
computed for the 5 % central A+A collisions. 

The pQCD calculation of minijet production is a momentum space calculation. In 
order to define the initial densities, a connection between the momentum of the minijet 
and its space-time formation point is needed. At collider energies the nuclei are strongly 
contracted, 2i?^/7 cm <C 1 fm, and the typical partons at saturation have a longitudinal 
spread of 2f(x\fs) ~ l/p S at 1 fm. We therefore consider the collision region as a 
point in the longitudinal direction and assume that the rapidity of the minijet coincides 
with the space-time rapidity of the formation point, y — r\ — (1/2) ln[(t + z)/(t — z)}. 
The formation (proper) time we take to be the inverse of the saturation scale, tq = 
1/Psat- Thus the minijet matter forms along the hyperbola t = J z 2 + Tq with initial 
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y/s/A [GeVJ 


130 


200 


5500 


Centrality 


% 


5 % 


r\ Cr/ 

% 


5 % 


-i r\ Cr/ 

10 % 


-i i*- Cr/ 

15 % 


r\ Cr/ 

% 


5 % 


A, 


197 


181 


197 


181 


166 


153 


208 


193 


T AcsAcB (0) [1/mb] 


29.4 


26.1 


29.4 


26.1 


23.1 


20.6 


31.7 


28.5 


■''part 


377 


346 


379 


347 


318 


293 


407 


376 


Psat [GeV] 


1.08 


1.06 


1.16 


1.15 


1.14 


1.12 


2.03 


2.01 


t [fm/c] 


0.18 


0.19 


0.17 


0.17 


0.17 


0.18 


0.10 


0.10 


a(E T ) [mbGeV] 


65.1 


67.0 


83.6 


85.9 


88.3 


90.7 


468 


479 


a(N 9 ) [mb] 


35.9 


37.3 


42.9 


44.6 


46.4 


48.2 


135 


140 


£,a<jV*> [mb] 


4.50 


4.66 


4.43 


4.58 


4.75 


4.92 


6.46 


6.66 


Eq<r(N q ) [mb] 


2.65 


2.76 


2.87 


2.98 


3.10 


3.21 


6.14 


6.34 


o-(iV) [mb] 


43.1 


44.7 


50.2 


52.2 


54.3 


56.3 


148 


153 


a(N B _ B ) [mb] 


0.616 


0.635 


0.520 


0.536 


0.553 


0.570 


0.106 


0.109 



Table 1: The pQCD quantities computed from Eq. Q with po = p sa t in central A+ A and in 
central A e g+A e g collisions at RHIC and LHC corresponding to various centrality cuts and 
Ay = 1. The saturation scales from Eq. (^Q), the corresponding formation times, and also the 
number of participants and the overlap function at b = are shown. 



longitudinal flow velocity v z (tq) = z/t. The partons in the rapidity interval Ay thus 
occupy a volume AV(r ) = A T Az = irR\T Ay at proper time r . 

The average initial energy density, number density and, as a new element |4TH Iq"2]. the 
net-baryon number density of the produced parton matter in central A+A collisions 
are then obtained as 

(e) = A y^ E T>AA {0, Paat ,Ay,y/I) (9) 
(n) = A ^ ^A(0,Psat,Ay,v^) (10) 
M = ^-jN^f(0,p sat ,Ay,^). (11) 

For nuclei with realistic transverse profiles, however, the saturation scale and thus 
also the formation time may vary with the transverse coordinate r = (x, y) jHHl EH] ■ 
We avoid such further modeling here and, following Ref. we consider a constant 
formation time ro = l/p sa .t and extract the transverse profile for the initial densities 
by differentiating the overlap function T A b with respect to d 2 r. The nucleon-nucleon 
luminosity for a transverse area element d 2 r is T A (h — r)Te(r) and the volume element 
dV = dzd 2 r = rAyd 2 r. These lead to the r dependent energy density and net-baryon 
density at (r = l/p sat ,r] = 0) in central (b = 0) A+A collisions, 

£ (r ,r) = J^ r= T A (r)T A (r f {ET) £>""' (12) 
T dr]d z r T Ay 
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dN AA B a(N B • B )Ay, Psat / 1Q \ 

n B (r ,r) = — — — = T a (t)T a (t) . 13 

T ar]d z r T Ay 

In Fig. we show the initial energy density in 5 % most central Au+Au collisions 
at v/i = 200 AGeV at RHIC and Pb+Pb collisions at ^/s = 5500 AGeV at the LHC. 
Notice the very large densities in the very center of the system. It should, however, be 
emphasized that since the volume near r ~ is small, the average energy densities are 
clearly smaller than those in the center, see Table 1 of Ref. [T2j. Another interesting 
point is the large difference between RHIC and LHC. This is on one hand due to the 
enhanced gluon production at higher energies (see Table HJ) and on the other due to the 
shorter formation time (smaller initial QGP volumes) in the LHC case (see Eq. Q12|)). 
Since the initial transverse size is nearly the same in both cases, the radial pressure 
gradients are larger at the LHC leading to stronger flow, as discussed in [60 . 
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Figure 1: Transverse profile of the 
initial energy density in an average 
A+A collision in 5 % most central colli- 
sions at the LHC (solid line) and RHIC 
(dashed lines) . Formation times and ef- 
fective nuclei are indicated, see Tabled 



The net-baryon number in central rapidity unit, Eq. (jHJ), can now be obtained from 
Tabled For RHIC, we get N B _ B = 16.6 and 14.0 in 5 % most central Au+Au collisions 
at ^snn = 130 and 200 GeV, respectively. For the LHC, we predict N B __§ = 3.11 in 
5 % most central Pb+Pb collisions at ^snn = 5500 GeV. 



3 From dense initial state to decoupled hadrons 
3.1 The hydrodynamical framework 

For the initial state of matter, the pQCD+saturation model provides the initial energy 
density e and net-baryon number density n B over a transverse profile at the formation 
time l/p sa t of the system. We make a bold assumption of fast thermalization by 
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taking the thermalization time scale to be the same as the formation time, tq = l/p sa t- 
Some arguments can be given in favour of this assumption. First, as discussed in 
[T7[ IdT] , from the point of view of (e) and (n) , the saturated (fully gluonic) system 
indeed looks thermal. This suggests that only elastic collisions causing energy and 
momentum transfer are needed for thermalization, number-changing reactions are not 
essential. Second, the formation time of more energetic (mini)jets is shorter than the 
time scale r associated with the saturation scale, so some secondary interactions have 
taken place already before the overall formation time. Third, recent studies of plasma 
instabilities lend support to fast thermalization [U21 EHj instead of the slower bottom- 
up thermalization jHlj. We also note that we assume full chemical equilibrium and 
follow the evolution of net-baryon number, which as a conserved quantity should not 
be sensitive on the details of chemical processes. We do not discuss the dynamics of 
chemical equilibration of the QGP here. Recent studies on this can be found in Ref. 

3.1.1 Simplifications: longitudinal scaling flow and azimuthal symmetry 

As discussed in Sec. 12.21 in our minijet description the parton system forms with an 
initial longitudinal flow velocity tanh# = v z = z/t. Assuming that the expansion main- 
tains this form, the original (3+l)-dimensional problem reduces to (2+l)-dimensional 
evolution equations for transverse expansion. For central collisions the final state is 
azimuthally symmetric leaving only the (l+l)-dimensional equations for radial flow to 
be solved. In this study we do not consider non-central collisions which exhibit the 
interesting phenomena of azimuthally asymmetric flow. 

Obviously boost invariance will be broken in real collisions with finite collision en- 
ergy. However, the assumption of boost invariant flow is still very good in the central 
rapidity region, as will be argued next. The effects of longitudinal pressure gradients 
on longitudinal velocity were studied in [37] , where the initial densities were computed 
from the pQCD minijet rapidity distributions and the initial longitudinal flow velocity 
was assumed to satisfy the scaling form. Since the radial expansion was neglected, the 
study I27J overestimates the change of the flow because cooling is slower and system 
lives longer than with radial flow. The obtained flow deviated only marginally (by a 
few per cent) from the scaling flow at the end of the evolution in the central rapidity 
region. Similar observation was made in a (2+l)-dimensional analysis at RHIC energy 

m 

As the scaling flow does not change much, the transfer of entropy, deposited ini- 
tially to the central rapidity unit, to larger rapidities due to the development of the 
longitudinal flow is negligible. Also, the dominant part of the final (thermal) spectra 
of hadrons comes from the part of the decoupling surface where \y — 9\ < 1, because 
only a small fraction of particles in thermal bath have large enough energy to end up 
in central rapidity over a rapidity gap larger than one. All in all this means that the 
fragmentation regions, where the assumption of boost invariance breaks down, do not 
affect significantly the results obtained using scaling flow at the central rapidity. 
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The question of a rapidity plateau in the final rapidity density dN/dy of hadrons 
is discussed in more details by Hirano et al. [SHI- They show that the plateau would 
appear in the measured spectra only if the initial densities are flat over a considerable 
range in space-time rapidity. At RHIC the measured pion spectra are Gaussian (HE] 
and it remains to be seen if a plateau develops at LHC energy. In any case, as argued 
above, a measured Gaussian rapidity density does not indicate that the description 
of transverse flow in the central rapidity region based on the assumption of scaling 
longitudinal velocity is an invalid approximation. 

We note that the approximations mentioned above can be relaxed in different ways. 
There has been a lot of discussion on elliptic flow, see e.g. a review jHZ], in which case 
the approximation of boost invariant flow is retained but the azimuthal symmetry is 
relaxed. Also, as mentioned above, a genuinely 3+1 dimensional hydrodynamical code 
for collider energies has been successfully developed jHHl EM ED] • 

3.1.2 Initial transverse flow? 

To complete the initial conditions for the hydrodynamic calculation the initial collective 
transverse motion, f r (r , r), must be specified. It can be argued that the change in the 
number of primary collisions, when going from the center to the edges in the transverse 
plane, could lead to non-zero initial transverse flow. However, hydrodynamical studies 
at the SPS energy show that only a very weak initial flow can be tolerated without 
getting into disagreement with the slopes of transverse hadron spectra [71]. Since 
there is no obvious reason for a strong transverse initial flow, we take for simplicity 
v r {r ,r) = 0. 

In Ref. [72] it was found that especially the spectra of protons, at RHIC, are better 
reproduced if some initial velocity is introduced. However, they use tq = 0.6 fm/c 
as the thermalization time and argue that the initial velocity arises from pre-thermal 
interactions. The velocity they assume is similar in magnitude to that created in 
hydrodynamical evolution up to r = 0.6 fm/c if r = l/p sa t ~ 0.2 fm/c is assumed, as 
we do. Whether or not the system is fully thermalized at r = l/p sa ,t, this early starting 
of hydrodynamical evolution can effectively describe the initial buildup of transverse 
flow as a result of the collisions which lead to thermalization. 

3.1.3 Equation of State 

Another key input in the hydrodynamic studies is the Equation of State, P = P(e, n#, ris), 
specifying pressure as a function of energy density, net-baryon number density and 
net strangeness density. Since strong interactions conserve strangeness, the total 
strangeness in the expanding matter of a heavy ion collision is zero. In principle 
non-zero strangeness spots could be generated in the matter through fluctuations but 
we assume that the net strangeness vanishes locally. This condition is imposed on the 
EoS and we write simply P(e,n B ) for P(e,riB,n s = 0). With the initial conditions 
considered here, the effects of n# on pressure are negligible, which we have checked, 
and for the hydrodynamical simulations we use P(e,riB = 0) for simplicity. However, 



11 



when calculating the final hadron spectra the value of tib and the condition ns = are 
essential in determining the baryon and strangeness chemical potentials which enter the 
thermal distributions of hadrons (see Eqs. (fTtjj) and (JT7J) ahead). When hydrodynamic 
simulations are converted to hadron spectra, the EoS in full (e, ne)-plane is used. 

In principle the EoS can be obtained from lattice QCD. A full EoS with finite quark 
masses and finite baryochemical potential fis is, however, not yet available. It is 
expected that for small values of /i# relevant at central rapidity region at collider 
energies the transition from hadron gas to parton matter is a rapid crossover j7Hj. For 
the recent progress on lattice calculations, see e.g. [73JEI]. From the point of view 
of radial hydrodynamic flow the difference between a weak first order transition and a 
rapid cross over is qualitatively not very significant as long as the EoSs are similar away 
from the transition region. Since we would like to be able to include the dependence of 
spectra on quantum numbers through the chemical potentials, we use in calculations 
an EoS with ideal quark-gluon plasma (QGP) in the high temperature phase and a 
hadron resonance gas (HRG) below the transition. For these phases it is easy to include 
the chemical potentials into the EoS. 

A first order transition is implemented by introducing a bag constant B into the QGP 
phase and connecting the two phases via the Maxwell construction. The inclusion of 
resonances in the hadron phase mimics the effects of both attractive and repulsive 
interactions between hadrons reasonably well [TH]. We include all hadron and hadron 
resonance states up to m = 2 GeV. A detailed account of constructing such an EoS 
can be found e.g. in Ref. [221 • With Nj = 3 we set the bag constant B to be B 1 / 4 = 
239 MeV, giving T c = 165 MeV for the transition temperature. 

Another way to construct an EoS would be to parametrize the n# = lattice results 
above the QCD phase transition [77] and join it smoothly with HRG EoS at T c . We have 
checked that hadron px spectra obtained with such an EoS do not differ significantly 
from those obtained with the bag model QGP HRG EoS 2 . In this paper, we shall show 
results with the QGP-HRG EoS only. 

3.1.4 Hydrodynamic equations 

With the longitudinal scaling flow and the assumed azimuthal symmetry, the transverse 
velocity and thermodynamic densities become functions of proper time r and transverse 
coordinate r only: v r = v r (r, r), e = e(r, r), n B = n B (r,r), P = P(r, r), .... The 
hydrodynamic equations are 

d,T^(x) = and d fl f B (x)=0, (14) 

where 

T" u = (e + P)u»u v - Pg^ and f B = n B u^ 

are the stress-energy tensor and the net-baryon current. With the above assumptions 
equations (|14j) simplify to 1+1 dimensional forms f01[7Jl]. Given an EoS in the form 
P = P(e, n B = 0) with n s = as discussed above, they can be easily solved numerically. 

2 For a recent study on the EoS effects on elliptic flow, see |75) . 
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3.1.5 Decoupling and decays 



In hydro dynamic description the final state hadron spectra are obtained by folding the 
thermal motion with the fluid motion on a decoupling surface a(r,r), which we define 
from the condition T(r, r) = Td ec - Given the velocity field and the thermodynamic 
densities, we calculate the thermal spectra of hadrons h by using the Cooper and Frye 
prescription [80 j 



dN h dN h 



d 3 p/E dydp\ 

9h_ 
2ir 



7T / da ll (x)p tl f h (x,p;T(x),iJ, h (x)) (15) 

J u 

CO . 

5](±l) n+1 / rTe ^(r,r)/T(r,r) ^ 
n=l 

[ -p T h{n^ r v r Y)^o(^7r^f) dr + m£/ (n7 r tv^)ifi(n7 r ^) dr] , 
where fXh = Bh/J-B + Sh^s is the chemical potential, gh the spin degeneracy factor and 



rrij, = \Jm\ + p\ the transverse mass of the hadron of type h. The last form, Eq. (|16|). 
holds for the boost invariant cylindrically symmetric radial flow. The upper (lower) 
sign is for mesons (baryons) and 7 r = (1 — v 2 )^ 1 ^ 2 . Also, the condition /^(r, r) < 
must always be satisfied for mesons, and, for the expansion in Eq. (JTHJ) to hold, also 
for baryons. 3 This is the case here. 

Chemical potentials for baryon number and strangeness, //g and are obtained 
on the decoupling surface by expressing the calculated net-baryon density and the 
strangeness neutrality in terms of free-particle densities: 

n B = J2 B hn h (T,Hh) and = ^ S h n h (T, (17) 

h h 

where Bh is the baryon number and Sh the strangeness of a hadron h. 

The preferred value Td ec = 150 MeV is extracted on the basis of RHIC data 
but we will also show the sensitivity of our results to T^ ec . The thermal spectra are 
calculated for all hadrons and resonances included in the hadronic EoS. 

To obtain final spectra which can be compared with the experimentally observable 
ones, we consider all strong and electromagnetic two- and three body decays of hadron 
resonances, with branching ratios obtained from [HI]. Decay kinematics is discussed 
in details e.g. in jH2J. Also feed-down effects from weak decays of hyperons will be 
discussed. 



3.2 Properties and effects of collective motion 
3.2.1 Longitudinal flow and transfer E T 

The effects of collective motion arise from the thrust due to collisions among produced 
particles, described by pressure in the hydrodynamic treatment. The most important 

3 For Hh(i~, r) > rrih, regions with Hh(i~, r) bigger or smaller than m\ must be treated separately. 
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overall effect of the rapid longitudinal expansion is to transfer a large fraction of (trans- 
verse) energy from central region to longitudinal motion, see Ref. jSZj for a discussion. 
The measured final transverse energy at RHIC is roughly a third of what is predicted 
by pQCD+saturation, and the same is expected to happen at the LHC as well [3]. At 
the same time, however, the number of calculated minijets is quite close to the final 
multiplicity. In the hydrodynamic approach the approximate conservation of the num- 
ber of particles is a consequence of the (approximate) conservation of entropy. The 
initial energy density, or equivalently, since instant thermalization is assumed, the ini- 
tial entropy density obtained from pQCD+saturation thus essentially determines the 
final multiplicity |17j . 

Longitudinal expansion alone would cool and dilute the matter to free hadrons. 
However, with the very large initial densities created at collider energies, the time 
needed for this, the longitudinal time scale, is so long that effects in the transverse 
direction cannot be neglected. From the point of view of the transverse energy at 
central rapidities we have two opposing effects: In the longitudinal expansion the large 
initial thermal energy is transferred to longitudinal motion and the transverse energy 
is reduced. On the other hand, the build-up of transverse flow limits this transfer and 
part of the initial thermal energy is converted to the transverse collective motion and 
thus to the measured transverse energy. Even though the fraction of energy transfer 
into transverse collective motion is much less than into longitudinal, it produces an 
effect which is seen as an increase in the average transverse momentum (see e.g. [83]) 
or in the effective temperature, measured as the inverse of the slope of the spectra. 

3.2.2 Properties of transverse flow 

The main overall effect of transverse flow in the hydrodynamic evolution is a more 
rapid cooling, leading to clearly shorter lifetime of the thermal system than with the 
longitudinal expansion alone. This is illustrated in Figs. [2HS1 which show our results 
for the hydrodynamic spacetime evolution in the framework discussed above. The 
figures E] and [3] are for RHIC and the next two for the LHC. The figures [2] and [3] 
show the flow lines in the (r, r) plane, and the figures [3] and [3] the transverse velocity 
contours. The phase boundaries QGP-mixed phase (MP) and MP HRG are shown by 
the thick curves. Also drawn are two isotherms corresponding to T& ec = 150 MeV and 
T dec = 120 MeV. 

It is interesting to compare the QGP lifetime at r = with that in a system subject 
to longitudinal scaling flow only. Such a system would cool down from €q to the critical 
energy density of QGP, e c = 1.93 GeVfm" 3 at r c = r (e /e c ) 3 / 4 . For the maximum 
values from the Fig. HJ e « 220 GeVfrn" 3 at RHIC and 2260 GeVfm" 3 at the LHC, 
and for the initial times from Tabled we get r c = 6.0 fm/c at RHIC and r c = 20 fm/c 
at the LHC. These are to be contrasted with r c (r = 0) ~ 5 fm/c at RHIC and 7.5 fm/c 
at the LHC in Figs. EJ . .0 The effects of transverse expansion thus become important 
already at the QGP phase at the LHC. 

The evolution of the transverse flow is clearly seen in these figures. The flow lines 
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Figure 2: The phase boundaries and the Figure 3: Contours of transverse velocity v r 

two decoupling curves studied, Td ec = 150 in the (r, r)-plane in 5% most central Au+Au 

and 120 MeV, with selected flowlines in the collisions at y / J/\F/v = 200 GeV. Reading from 

plane of proper time r and transverse vari- left to right, v r = 0.1, 0.2, . . . ,0.9. 
able r in 5% most central Au+Au collisions 
at ^/s]^n = 200 GeV. 
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Figure 4: As Fig. but for 5% most central Figure 5: As Fig. El but for 5% most central 
Pb+Pb collision at ^/sWn = 5500 GeV. Pb+Pb collision at ^/sWn = 5500 GeV. 
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r = r(t) (at z = 0, r = i) in Figs. 121 and H] depict how a fluid element in a given 
differential interval in transverse variable at r moves radially outwards. There is no 
entropy-flux across the flow lines so we get an idea of where the initial multiplicity 
goes (once the transverse volume is properly accounted for). The slope of the flow line, 
dr(t)/dt, gives the local radial flow velocity v r (r(t),t). Thus, the bending of the flow 
lines seen in the QGP and HRG phases indicates acceleration of the fluid. The straight 
flow lines observed in the MP in turn indicate a constant pressure and no acceleration. 

The corresponding constant velocity contours at RHIC and LHC, which, reading 
from left to right are v r = 0.1, 0.2, . . . , 0.9, are shown in Figs. El and The velocity 
contours illustrate how the deflagration wave penetrates into the fireball as the radial 
pressure gradient pushes the matter to collective motion. A kink at the boundary of 
QGP and mixed phase (MP) results from disappearance of pressure gradients in the 
MP. As there is no acceleration in the MP, the matter continues to flow with radial 
velocity obtained already in QGP. The reappearance of pressure gradients, when matter 
enters from MP to hadron gas phase, is seen as a kink in the velocity contours at the 
phase boundary. 

The dense bundle of velocity contours starting around r > 10 fm results from the 
treatment of the nuclear surface. Even though we take the densities smoothly to zero, 
the expansion is essentially to vacuum and the velocities approach the velocity of light. 
This region has no consequences for any physical observables, since it is outside the 
decoupling surface in a region with vanishingly small densities. 

Since the larger energy density leads to larger pressure gradients at LHC than at 
RHIC, the transverse flow at LHC grows also much stronger than at RHIC. This is seen 
both in the larger radial extent of the decoupling curve and the larger values reached 
in transverse velocity at decoupling. We next discuss how this affects qualitatively the 
transverse spectra of hadrons. 

3.2.3 Effects of transverse flow on hadron spectra 

It is well known that transverse flow broadens the ttit spectra of hadrons and, since 
heavier particles obtain a larger momentum increase from the collective flow velocity, 
the effective temperature T e g, defined as the inverse of the slope of the logarithm 
of a spectrum, and the average transverse momentum 



increase with particle mass. Before discussing the actual results, we demonstrate the 
effects of flow in a simplified system of homogeneous cylinder with radius R, which 
decouples at r = t/ = const. In this case, if Hh = 0, the Cooper-Frye formula (|T5|) 



(Pt) 



J dpx Pt dN h / dydp T 




J dp T dN h I dydp T 
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reduces to 4 



" A fr f nvr E(±l) n+1 f* ' dr rl {n^v^K^n^), (19) 
2tt n=1 Jo T T 



dydm 



2 



where the index h has been suppressed. Integration over r is trivial if v r = const. In 
particular, if no transverse flow is present, v r = 0, the equation (|T9"|) simplifies to 

dN 9 n2 ^/MW+lr//™^™T/™ 9 



^^Et^r^^J^-^^V^re-^, ,20, 

where also the asymptotic behaviour at large mx/T is indicated. 

The upper set of curves in Figgis obtained using Eq. (|20|1. i.e. there is no transverse 
flow. In that case, spectra of all particles, here three spin-0 bosons with masses m = 
140, 494 and 938 MeV, are the same except that their m^-spectra start from particles 
mass (different curves are made visible by multiplying the m = 494 and 938 MeV 
curves with 1.2 and 0.8, respectively). This is known as the m^-scaling. The overall 
normalization of the curves is left free, as we consider now the slopes only. 

Non-zero transverse velocity breaks the m^-scaling at pr < m. This is illustrated by 
the lower set of curves in Fig. fusing constant transverse flow v r = 0.6. The asymptotic 
behaviour of Eq. (fT9~|) at itlt/T ^> 1 shows an exponential factor exp{— my/T e ff } with 



T e s = ^jz^T. Thus the scaling of the slopes is recovered asymptotically. In the 
example of Fig. Hwith v r = 0.6, T e R = 2T. This explains why the slopes of the spectra 
are nearly the same in the cases T = 220 MeV, v r = and T = 110 MeV, v r = 0.6. 

This simple study very clearly demonstrates the role of the transverse flow in the 
properties of transverse spectra of hadrons, and the necessity to understand the dy- 
namical evolution of the transverse flow in detail for proper interpretation of measured 
slopes. 

Let us turn to discuss the flow in more realistic hydrodynamical simulations, pre- 
sented in Figs. EHSl There is an interesting interdependence between the EoS, space- 
time evolution of the flow and the slopes of thermal hadron spectra. As long as a 
phase transition is included in the EoS, the strength of final flow on the decoupling 
surface, of fixed decoupling temperature, is not very sensitive to the number of hadron 
resonance states included in the hadronic phase, while the actual spacetime evolution 
of the system may vary quite significantly. This can be understood as follows: Re- 
duction of the hadronic states in the EoS leads to a stronger phase transition between 
the QGP and HRG, i.e. latent heat as well as entropy ratio between the two phases 
increase. This makes the deflagration from mixed phase to HRG stronger. It increases 
also the pressure gradients in HRG since the reduction of hadronic states leads to a 
larger sound velocity c s and the pressure is given by p = c 2 s e. On the other hand, 
as discussed in [52] . the energy release from heavy resonance states significantly slows 
down the decrease of temperature during the evolution. Thus, when we decrease the 

4 This form, Eq (|19fl . was used in a blast-wave fit |84| . where v r — (r/i? J 4)u max was assumed and 
where the parameters T and i> max were fitted to reproduce the shape of the final hadron spectra. 
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Figure 6: An example of how trans- 
verse flow affects the slopes of hadron 
spectra in a simplified system which 
decouples at a constant proper time 
and a constant radial flow velocity v r . 
Two sets of spectra, computed from 
Eq. (|19j) but leaving the overall normal- 
ization arbitrary, are shown for bosons 
of masses m = 140 (solid lines), 494 
(dotted) and 938 MeV (dashed). The 
upper set of curves, where v r = and 
temperature T = 220 MeV, shows rriT 
scaling (for better visibility, the curves 
for 494 and 938 MeV are multiplied 
by 1.2 and 0.8, respectively). In the 
lower set v r = 0.6 and T = 110 MeV. 
At small rriT, the scaling is broken 
due to a non-zero v r , but the spectral 
slopes are similar to the upper set at 
rriT 3> T, m. 

number of hadron resonance states in HRG, it results in faster development of flow and 
drop in temperature, leading to a shorter lifetime of the matter. For larger number of 
hadron resonance states the flow grows more slowly but the energy release from heavy 
resonances leads to longer time for the flow to build up. In practice these effects com- 
pensate each other and as a result, the slopes of the spectra are relatively independent 
of the hadron content of the EoS. 

There is, however, a relatively strong dependence of the effective temperature, the 
inverse of spectral slope, on the decoupling temperature [22] once the hadron content 
of the EoS is fixed. In Sec. El we shall illustrate the variation of the hydrodynamical 
results in the range of decoupling temperatures T& ec = 150 . . . 120 MeV. In particular, 
with the given pQCD+saturation initial state and with the relatively steep initial 
transverse profile, we shall show that the RHIC spectra favour the higher T^ ec both at 
y/s = 130 GeV and 200 GeV. 

4 Hadrons from fragmentation of high-py partons 

4.1 pQCD framework for computing high-pr spectra. 

Inclusive hadron spectra in high energy hadronic or nuclear collisions are computable in 
terms of collinearly factorized partonic hard parts (cross sections in LO) &ij^f + k, par- 
ton distribution functions (PDFs) fa/ a, fj/B and fragmentation functions (FFs) Df^, 
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schematically as 



da AB-, h+ x = j2 f i/A g f. /B g g D/^, (21) 

where the fragmentation function, Df_>h, gives the average multiplicity of h hadrons 
produced from a parton /. The NLO pQCD computation [HEl IBE1, which requires 
NLO PDFs and FFs as well, provides a successful description of the high-p^ spectra 
in p+p collisions at collider energies. Also the leading-order considerations are useful 
and have been shown to agree quite well with the data in the large-py region, once 
the overall normalization, the K factors, are extracted from the data ^B]- We shall 
here make use of the latter: by applying the energy- dependent K factors of Ref. jTHj 
and including in the analysis also the recent p+p data from RHIC jS], we determine 
the K factors at \fs = 130, 200 GeV and also, through extrapolation, at the LHC 
energy y/s = 5500 GeV. Armed with these and by including nuclear effects in the 
PDFs according to the EKS98 parametrization we compute the baseline pQCD 
spectra of high-p^ hadrons for nearly central A+A collisions at RHIC and LHC. We do 
not make further modeling to account for the Cronin effect, or any intrinsic effects 
in the PDFs or in the FFs [EH EH El ■ 

As presented in detail e.g. in Ref. ^H], inclusive spectrum of the high-p^ hadrons 
in A+A collisions at an impact parameter b in the collinearly factorized LO pQCD 
formalism becomes 



dNhAh) r v , dz ^ . 2 .dN{ A {b) 

= K(VS, Q, mp) " J(mr, v) E F D„(z, A)-^f 



[22) 



where f = g,u,d, ... labels the different parton types and h those of hadrons. Transverse 
momentum and rapidity of the produced parton are qT > qro and Vf, while and y 
are the transverse mass and rapidity of the hadron. The hadron of fractional energy 
z = Eh/Ef is assumed to form collinearly with its mother parton. Thus the kinematic 
variables are related as 

q T = ^J( mT ,y), J(mr,y) = (l- 2 m \ 2 )~ 1/2 (23) 
z m T cosh y 

yj = arsinh( — -sinhy), and — ^ cosh y < z < min[l, ^-J(jriT, y)]. 
Pt qTo 



We shall here apply the KKP fragmentation functions [90J with a scale choice mf = Pt- 
The inclusive spectra for the production of parton type / is obtained at LO from Eq. (JSJ) 
by integration 



dN f AA (b) da f AA _ T , h ) f dydy y do_ AB ^ !+s 

dq T dy f AA dq T dy f AA J 1 2 ' ^ dq T dy ± dy 2 

1 



[SkfS(yf-yx)+S lf 5(y f -y 2 )\——, (24) 

J- + Okl 
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where Taa is the nuclear overlap function from Eq. (JJJ). We define the nuclear PDFs 
by including the x- and Q-dependent EKS98 nuclear modifications jjjj on top of the 
the CTEQ5 (§3 PDFs of the free proton. Isospin effects in the nPDFs are accounted 
for in the same manner as in Ref. jUj. The factorization/renormalization scale is again 
chosen as Q = qr- 

The factor K(s/s, Q, /xp) controlling the overall normalization of the pQCD fragmen- 
tation spectrum in Eq. (|22|). is here a phenomenological parameter which accounts for 
all higher-order contributions in the partonic hard parts, PDFs and FFs (as well as 
for any scale dependence arising from a truncation of the perturbation series). This K 
factor can be determined on the basis of the p+p and p+p data. It should be empha- 
sized that K depends not only on y/s, but also the choices for the scales Q and /xp, and 
on the PDF and FF sets used. The K factors for the parameters specified above have 
been obtained in Ref. ^BJ for various cms energies by using the data from the AFS 
|H2] . UA1 inn] , UA1 MIMI |nij and CDF jHS] experiments. These results are shown in 
Fig. [7| As new input, we now include also the K factor corresponding to the recent 
STAR charged particle data from p+p collisions at y/s = 200 GeV at RHIC 6J. The 
inner error bars shown correspond to the statistical errors of our i^-factor analysis JH] > 
for the outer error bars these have been added in quadrature with the systematic errors 
of the data 5 . As discussed in detail in JTH] , the best agreement between the shapes of 
the computed and measured spectra is obtained by construction at pt ^ 4 . . . 5 GeV 
but the overall agreement is surprisingly good over a wider p? range. 

Based on the systematic decrease of K(y/s) in Fig. [7J we extrapolate our LO pQCD 
reference spectra to the LHC by making the following two simple parametrizations to 
the obtained i^-factors, 



where y/s is in GeV's. We get (a, b) = (3.4131; -0.44569) for the fit I with x 2 = 43.2, 
and (a,b,c) = (5.6206; -1.2139; 0.065615) with X 2 = 40.6 for the fit II. Both fits are 
also shown in Fig. Since systematic errors are not given for all cases, only statistical 
errors (see |18j ) are considered in these fits. At RHIC energies the K factors are not 
sensitive to the assumed form of parametrization, see Fig. [7| For the pQCD spectra at 
the LHC, we shall give an uncertainty band corresponding to the different i\T-factors 
obtained from these fits. We note that a similar fit procedure was introduced also in 
but modeling in also a Gaussian smearing from intrinsic-fcr effects. 

4.2 Energy losses for hard partons 

To study the interrelation of the hydrodynamic and pQCD hadron spectra at RHIC and 
at the LHC more meaningfully, we also consider effects induced by the partonic energy 
losses to the high-p^ pQCD fragmentation spectra. First we note that the nuclear 

5 For CDF and UA1 MIMI we assume a 15 % systematic error. 
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a + b In a/s 
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Figure 7: The K factors obtained from fitting the computed leading-order high-p^ hadron 
spectra to the p+p and p+p data. The fit I to lniT, which is linear in lni/i, see Eq. ()26jl . 
is shown by the solid line and the quadratic fit II by the dashed curve. The parameter 
values and error bar specifications can be found in the text. 



effects in the PDFs, which are constrained to reproduce the nuclear DIS data through 
global DGLAP fits, are not responsible for the factor ~ 5 suppression discovered in 
the hadron spectra in central Au+Au collisions at RHIC at pr^5 GeV. In fact, as 
shown in JH], the nuclear modifications slightly (by 0...15%) enhance the spectra below 
Pt ~ 13 GeV relative to p+p. The energy losses we consider within the framework 
given in Refs. fHl 120]; by folding quenching weights with parton production cross 
sections and FFs, schematically written as 

da AA^ h+ x = ^ f , /A ^ fj/A 3 £.._^ +fc Pf ( AEi L; q ) g, Df ^ (27) 

ijk 

The quenching weights, Pf(AE, L, q), are generalized probabilities for a hard parton / 
to lose an amount AE of its original energy through radiation induced by the medium 
over a path length L. The transport coefficient q corresponds to the average transverse 
momentum squared per unit pathlength gained by the hard parton in traversing the 
matter. The transport coefficient reflects the average number and energy densities of 
the system as q oc n ~ e 3 / 4 [96J, and it controls the energy loss. 

We compute the pQCD spectra subject to energy losses according to the formulation 
given in sec. 2.2 of Ref. [TH]. We do not repeat the discussion here, but for illustrative 
purposes we note that when the energy losses are taken into account, the analogue of 
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Eq. (}2*2*j) can be written in terms of medium-modified fragmentation functions [HZj for 
central A+A collisions as 

Vgfi = K ■ A m T , *r*lffh r, ■ (28) 

where the hard parton's original transverse momentum is = J(rriT,y)pT/z, the 
distance of the hard parton's production point from the center of the transverse plane 
is r = |r|, and the azimuthal angle between the hard parton's original transverse 
momentum qx and the transverse vector r is <fi. The medium-dependent fragmentation 
functions are defined as 



I>KW.6j£) = / ^Pf(s,L(r,cf>),q)D f ^ h (^-y T ), (29) 



de 

-P f (e,L(r, <p),q)D f ^ h {— 

z 1 — 6 

where e = AE/Ei is the energy fraction lost by the hard parton / and the quenching 
weights are computed from those of Ref. |2()j . The length L traversed by the parton 
in the QGP now depends on the transverse location of production and on the angle 0. 
In determining the path lengths L, a homogeneous transverse disc profile of a radius 
Ra is assumed [T9"] . 

In what follows, we compute the spectra with normalization fixed by the i^-factors 
determined above. As discussed in [TU], the obtained spectra are not sensitive to the 
(average) lifetime of the plasma once it is of 4-5 fm/c or more, which is the case here, 
see Figs. El and |U In particular, in 5 % most central Au+Au collisions at ^snn = 200 
GeV we apply a large time-averaged transport coefficient q = c^oo — 10 GeV 2 /fm which 
was extracted in Ref. jTHj from the RHIC data in the region pt > 5 GeV. For the other 
energies and A e{i studied, we scale q according to 

q(A, v^7) = g 2 oo(A/181) - 383 (v^/200 GeV) ' 574 ( 30) 

on the basis of the average number density scaling of our QGP initial conditions, as 
discussed in (THj. 

The quenching weights applied here are computed in the eikonal high-energy ap- 
proximation where the lost energy AE can in principle be arbitrarily large. Thus, 
with finite kinematics it may happen that the probabilities are not always normalized 
properly to unity within the parton energy range available, 

pqt cosh y 

/ d(AE)P f (AE,q,L) < 1. (31) 

Although too large suppression can follow from this loss of probability, this computa- 
tion, however, gives us a lower limit of the hadron spectra. An estimate of the upper 
limit can be obtained by reweighting the quenching weights so that probability conser- 
vation is always enforced, see jHHl UH] • In our results for the high-p r spectra (Figs. |Hl 
El EH EU dni CHI CSl ahead) , we shall give an uncertainty band corresponding to these 
two limits. 
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5 Results 



We shall next compare the hydrodynamical and pQCD-based results with the mea- 
sured spectra in central and nearly central Au+Au collisions at RHIC and make an 
extrapolation to Pb+Pb at the LHC. At RHIC, the hydrodynamic spectra cover the 
small-pr region up to a few GeV and the pQCD spectra the region at Pt ^ 4 . . . 5 GeV. 
The pQCD calculation of the hadron spectra is expected to become less accurate at 
the lower end, below px < 4 GeV (see |TS|). and similarly the uncertainty in the tails 
of hydrodynamic results grows with px- However, as will be seen below, the differ- 
ence in the slopes of the spectra from the hydrodynamic and the pQCD calculations 
with energy losses in particular, is rather large. This should reduce the uncertainty in 
the crossing of two contributions at the LHC in particular, where the hydrodynamic 
spectra are shown to dominate over a larger px region than at RHIC. 

5.1 Comparison with RHIC data 

We first compare our results with the experimental data collected by STAR (HI 

unnj . phenix ihi isi nnn no2i nnsi nn^i , phobos [iniiinni and brahms nunnniiMi 

collaborations for the most central bins in Au+Au collisions at ^Jsnn — 130 and 200 
GeV. As explained in Sec. 12.21 the centrality selection is accounted for by considering 
central A e g+A e s collisions. 

5.1.1 Transverse momentum spectra 

First, we consider the px spectra of charged particles measured at central pseudo- 
rapidities rj. Fig. |H] shows the data taken in 5 % most central Au+Au collisions at 
v /5^7 = 130 GeV by STAR {ZJ and PHENIX gj collaborations. Similarly, Fig. H 
shows the spectra of charged particles at y / J/v r /v = 200 GeV measured by PHENIX jH] , 
STAR 0, PHOBOS HUJ and BRAHMS [H] Collaborations. The hydrodynamic results 
are computed accordingly [3], using two different decoupling temperatures, T dec = 150 
MeV (solid line) and 120 MeV (dashed line). In addition to the strong and electro- 
magnetic decays, we have here included also the feed down from the weak decays of 
hyperons, which affects the total number of charged pions, kaons and (anti)protons on 
a few-per-cent level (see Table HU ahead) . Note that our hydrodynamical results are 
for 5 % most central collisions while some data are taken in larger centrality bins, as 
indicated in the figures. 

Results in Figs. |H1 and El show that the normalization of the calculated hydrodynami- 
cal distributions is essentially independent of the decoupling temperature and since the 
slopes deviate only at larger transverse momenta, the total charged-particle multiplicity 
is not very sensitive to the decoupling temperature [3]. 

If the system continues to undergo hydrodynamical evolution at temperatures be- 
low 150 MeV, the collective flow will increase and this is seen as an enhancement of 
high p T particles. For y/s = 200 GeV, Fig. EE the calculation with T dcc = 120 MeV 
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Pt [GeV] 

Figure 8: Transverse momentum spectra of 
charged particles at 77 = (averaged over 
1 77 1 < 0.1) in 5 % most central Au+Au col- 
lisions at y/sJm = 130 GeV. Our hydrody- 
namic results are shown for T dec = 150 MeV 
(solid line) and T dcc = 120 MeV (dotted line). 
The pQCD fragmentation results are shown 
with energy losses (shaded band, see the text) 
and without (dashed line). The STAR data 
is plotted with the given total errors, and 
the PHENIX data jU] by adding the given sta- 
tistical and systematic errors in quadrature. 
Note the different centrality classes of the two 
data sets. 
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Pt [GeV] 

Figure 9: As Fig. |H]but for 5 % most central 
Au+Au collisions at ^/snn = 200 GeV. The 
data shown is taken by STAR 0, PHENIX 
[HI, PHOBOS QUI and BRAHMS Q3] in the 
centrality classes indicated. We have plotted 
the STAR and PHOBOS data with the given 
total error bars, the PHENIX data by adding 
the given statistical and systematic errors in 
quadrature and the BRAHMS data with the 
given statistical error bars. 



overshoots the experimental data already at px ~ 1 GeV, while a high decoupling 
temperature Td ec = 150 MeV describes the results better. It is also clearly seen that 
the hydrodynamic spectra cannot describe the RHIC data at larger transverse mo- 
menta, pt > 2 . . . 3 GeV. The same applicability region for the hydrodynamical spectra 
at RHIC has been suggested also by the analyses of azimuthal anisotropies, see e.g. 

Here it should be emphasized that we have not tuned the computation of the initial 
state (e.g. the possible unknown constant in the saturation criterion in Eq.QJ) to fit the 
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normalization to the data but kept it as in the original predictions fHIH]- Given the 
robustness and simplicity of the pQCD+saturation model, we obtain surprisingly good 
description for pion and kaon spectra. As pions dominate the multiplicity, our results 
for integrated observables, such as total multiplicity and transverse energy, are also 
close to the experimental results. Spectra of protons and heavier particles, especially 
at v /J/v r /v = 200 GeV, show details which are to be addressed in more detail. 

The dashed lines in Figs. |H] and |2] present the spectra of charged hadrons (charged 
pions, kaons and (anti)protons included) from fragmentation of partonic jets as pre- 
dicted by Eq. in the case of no energy losses. It is seen that this prediction is 
clearly above the data when grows beyond 4-5 GeV. The solid curves bordering the 
shaded area are the corresponding upper and lower limits of spectra computed with 
energy losses as explained in Sec. 14.21 In all these pQCD based curves, we apply the 
V^-dependent K factor according to fit II in Eq. (J2EJ). 

As observed earlier ^Hj and in these figures, the high-p^ spectrum is clearly sup- 
pressed relative to the baseline pQCD fragmentation spectrum without energy losses. 
When the energy losses are taken into account, the pQCD spectra agree with the data. 
Apart from the .ff -factor, this agreement is obtained by construction at ^/s^n — 200 
GeV since the value of the transport coefficient we use here was determined in ^2] by 
a fit to the data at px > 5 GeV. It is, however, quite interesting to see how the pQCD 
fragmentation+energy loss spectrum in turn gradually fails to describe the data at 
Pt ^ 4 . . . 5 GeV. Thus, in the intermediate-pr region, 2 ... 3 GeV < pr < 4 . . . 5 GeV 
at RHIC neither the hydrodynamic nor pQCD description is sufficient to describe the 
data alone but both are needed. We will return to this interesting crossover region and 
the independency of the two approaches in Sec. 15.31 

Next, we proceed to the pt spectra of identified hadrons at midrapidities. Fig. ITU1 
shows the PHENIX data collected for positive pions, kaons and protons in 5% most cen- 
tral Au+Au collisions for y = at = 130 GeV [M]. Similarly, in Fig. HJJ STAR 
1221, PHENIX HUH] and BRAHMS [MJ [M| data are shown at ^IJ^ = 200 GeV. 
Notice the scaling factors 10 and 1000 for kaons and protons, respectively. The hy- 
drodynamic spectra are computed with a common decoupling temperature, i.e. the 
thermal and the kinetic freeze-out take place simultaneously. The results with a high 
decoupling temperature, Td ec = 150 MeV, are shown by the solid curves, and those 
with a lower value, T^ ec = 120 MeV, by the dashed ones. It is seen that there is a 
substantial change in the normalization at small pt for protons and kaons, indicating 
also a change in their multiplicity dN/dy. For pions the change in the multiplicity 
with Tdec is small but the slopes change considerably for all three particles. Using a 
common high decoupling temperature Td ec = 150 MeV (solid lines) our results describe 
both the normalization and the slopes of the experimental spectra quite well whereas 
for Td ec = 120 MeV (dotted lines) the normalization of protons and kaons is too low 
and the pt dependence of all spectra is too shallow. 

The shaded bands in Figs. [TU] and [TT] show the spectra of pions, kaons (divided by 
10) and protons (divided by 1000) calculated from pQCD+ fragmentation+energy loss. 
For the pion spectra the hydrodynamic and the pQCD calculations complement nicely 
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Figure 10: Transverse momentum spectra 
of positive pions, positive kaons and pro- 
tons at y = in 5 % most central Au+Au 
collisions at ^/snn = 130 GeV. The solid 
lines show our hydrodynamic results with 
Tdec = 150 MeV and the dotted lines the 
results with T dec = 120 MeV. The shaded 
bands correspond to the pQCD fragmenta- 
tion results with energy losses. The PHENIX 
data |101j is plotted with the given total er- 
ror bars. Note the scaling factors 10 and 1000 
for kaons and protons, respectively. Both the 
hydrodynamic result and the PHENIX data 
contain the feed-down contributions from hy- 
perons. 



Figure 11: As Fig. E3 but at ^/snn = 
200 GeV. The PHENIX data QUI and the 
BRAHMS data pIMllMl are shown with sta- 
tistical errors and the STAR data jHH] by the 
given total error bars. The hydrodynamic 
calculation and the PHENIX data are with- 
out the hyperon feed-down contributions but 
the STAR and BRAHMS data contain the 
feed-down. 



each other: the former describes the data at low and the latter at high pr's- Even 
though tempting, it would be too naive to simply add the two contributions without 
further considerations. E.g. one would expect the thermalization assumption to fail 
in the large momentum tails of the thermal distribution modifying the region where 
the dominance of thermal part goes over to the dominance of the pQCD part. See the 
discussion in Sec. 15.31 However, for protons at pr ~ 3 GeV, the pQCD fragmentation 
contribution is clearly smaller than the hydrodynamic one and also clearly below the 
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measured proton spectra. We thus conclude that pQCD fragmentation contribution as 
considered here is not the origin of the intermediate-pT proton abundance. 

The PHENIX data |l(Jlj in Fig. El include the feed-down from weak decays of hy- 
perons 6 . Our hydrodynamic spectra are here shown with the feed-down as well. While 
for pions and kaons the hyperon feed-down is a small effect, for protons it corresponds 
to a ~ 30% increase. On the other hand, the hyperon feed-down is removed from the 
y/sjfij = 200 GeV PHENIX data PHH in Fi S- El but the STAR [H] and BRAHMS 
|106[ I107j data include feed-down contributions. Our hydrodynamic results are shown 
without feed-down. 

Our hydrodynamic antiproton and antilambda spectra without (dashed line) and 
with (solid line) the hyperon feed-down contributions are shown in Fig. ^] for 5 % 
most central Au+Au collisions at y/s^N = 130 GeV. The data in this figure are from 
STAR |100j (circles) and PHENIX |102j (triangles) in the corresponding centrality 
bin. Both data sets include the hyperon feed-down contributions and should thus be 
reproduced by the solid lines. Again, the agreement is surprisingly good, except at the 
largest px bins. 

Since neutral pions produce the most important background for direct and thermal 
photons 1081 11091 r7T] . an accurate description of neutral pion spectra would be par- 
ticularly important. In Fig. El we compare our tt° spectrum with that measured by 
PHENIX [inil in 10 % most central Au+Au collisions at y/s]^ = 200 GeV. Again, 
the hydrodynamic computation is shown for the low (dotted) and high (solid) Td ec , and 
the pQCD fragmentation+energy loss results, where we assume 7r° = (tt + + 7t~)/2, and 
the K factor from the charged-particle p+p(p) data (Fig. EJ), are shown by the shaded 
band. The conclusions here are the same as for charged pions: the high T^ cc describes 
the small-pr part well until px ~ 3 GeV, while the pQCD+energy loss spectrum re- 
produces the high-p-r part px > 5 GeV. In the intermediate px region, 3 < px < 5 GeV, 
both contributions are needed. 

Figure El shows experimental results for very low transverse momentum spectra of 
hadrons measured by the PHOBOS collaboration |105j . Also the spectra measured by 
the PHENIX collaboration 103J extending to larger px are shown. Both collaborations 
have removed feed-down contributions from their data. Our results for 7r + +7r~, K + +K~ 
and p+p are plotted as solid, dotted and dashed line, respectively. The measured 
spectra are quite well reproduced. On a log-scale the small-pr interval illuminates 
clearly that also spectrum of pions becomes flat when px < m n , which is also supported 
by the PHOBOS data. The spectra of kaons and protons we obtain, are slightly above 
the data. This will be seen in the multiplicities of these particles discussed in the next 
section. 

6 PHENIX has published a feed-down corrected spectrum of antiprotons in minimum bias events 
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Figure 12: Transverse momentum spectra 
of antiprotons and antilambdas at y = 
in 5 % most central Au+Au collisions at 
\/ s nn = 130 GeV. Our hydrodynamic re- 
sults are for Td ec = 150 MeV, the solid lines 
showing them with hyperon feed-down con- 
tributions and the dotted lines without. The 
PHENIX p [M] and A [TUH data and the 
STAR data 100 contain the feed-down con- 
tributions. In the PHENIX p data the error 
bars shown correspond to the given total er- 
rors, all others to the statistical errors. 



Figure 13: Neutral pion production in 10 % 
most central Au+Au collisions at ^/snn = 
200 GeV. Our hydrodynamic results are 
shown for Td cc = 150 MeV (solid line) and 
T dec = 120 MeV (dotted line). The pQCD 
fragmentation results are shown with energy 
losses (shaded band, see the text) and with- 
out (dashed line). The data is taken by 
PHENIX |104j and shown with the given to- 
tal error bars. 



5.1.2 Multiplicities and feed-down effects 

Experimental results for multiplicities, dN/dy, of n^, K^, p and p from PHENIX 
|103| and STAR [2| collaborations with ^/snn = 200 GeV are presented in Table El 
PHENIX has subtracted the feed-down from their data, but STAR data include it. 
Pion and kaon multiplicities in these experiments agree within the given errors. The 
difference between the PHENIX and STAR results for (anti)protons, however, may 
require a more detailed study of the feed-down effects. 

Table El summarizes our hydrodynamic results for integrated observables in 5% most 
central Au+Au collisions at RHIC energies \fs NN = 130 and 200 GeV. The results 
are quoted for the decoupling temperatures Td cc = 120 and 150 MeV. The total and 
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Figure 14: Transverse momentum spectra of charged pions (solid), charged kaons (dotted) 
and protons + antiprotons (dashed) at y = in 15 % most central Au+Au collisions at 
\/ s nn = 200 GeV in the ultralow-p^ region measured by PHOBOS |105j . Our hydrodynamic 
results are for T^ ec = 150 MeV. The PHOBOS data is shown with statistical and systematic 
error bars and the PHENIX data |103j by adding the statistical errors in quadrature. Note 
the different centrality classes in the two data sets. 



dN/dy 


7T+ 


7T 


K+ 


K" 


P 


P 


PHENIX 


286.4 ±24.2 


281.8 ±22.8 


48.9 ±6.3 


45.7 ±5.2 


18.4 ±2.6 


13.5 ± 1.8 


STAR 


322 ± 32 


327 ± 33 


51.3 ±7.7 


49.5 ± 7.4 


34.7 ±6.2 


26.7 ±4.0 



Table 2: Experimental results for dN/dy in 5 % most central Au+Au collisions with ^Jsnn = 
200 GeV from PHENIX [EH] and STAR [<S] Collaborations. The PHENIX results are 
corrected for the feed-down, in the STAR results feed-down is in the data. 

charged particle multiplicities and the total transverse energies at central rapidity are 
listed first. Below these, we show the total multiplicities for different hadron species. 
The last three rows of the table show our results for the antibaryon-to-baryon ratios 
for protons and lambdas, and the total net-baryon number. 

The first numbers correspond to results in which the hyperon feed-down contributions 
have not been considered, i.e. we include all the electromagnetic and strong decays but 
regard those hadrons stable which have only weak decay channels. In the results given 
in parentheses the hyperon feed-down contributions have been added. The feed-down 
from A, in particular, has been included to all the other yields while the feed-down for A 
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V s aw [LreVJ 


130 


200 


idee [MeV J 


120 


150 


120 


150 


dN/dy (tot) 


1086 


(1097) 


1141 (1178) 


1284 


(1297) 


1352 (1396) 


dN/dr)(tot) 


960 


(969) 


961 (988) 


1141 


(1152) 


1148 (1181) 


dN/dy (char) 


671 


(680) 


687 (716) 


794 


(804) 


815 (849) 


dN/dr] (char) 


597 


(604) 


587 (608) 


710 


(718) 


701 (726) 


dE T /dy(tot) [GeVJ 


696 


(696) 


709 (709) 


850 


(850) 


868 (868) 


drj^/arj [tot ) [*oreVJ 


625 


(625) 




769 


(769) 


(61 {(00) 


7 7V r7r+ / 7 


284 


(287) 


269 (279) 


337 


(340) 


319 (332) 


7 A T7T® / J 

aiv /ay 


312 


(316) 


298 (312) 


370 


(375) 


354 (370) 


dN n ~ /dy 


284 


(289) 


269 (282) 


337 


(342) 


319 (334) 


cW K /dy 


43.6 


(43.7) 


53.0 (53.1) 


51.3 


(51.3) 


62.4 (62.5) 


dN K ~/dy 


41.1 


(41.1) 


50.3 (50.4) 


49.1 


(49.1) 


60.1 (60.2) 


dN p /dy 


10.1 


(14.0) 


20.5 (30.2) 


10.4 


(14.5) 


23.0 (34.0) 


dN p /dy 

7 A T A / 7 

dN A /dy 
dN R /dy 


3.7 


(5.2) 


14.0 (21.1) 


5.0 


(7.1) 


17.5 (26.4) 


4.1 


(5.1) 


9.9 (12.9) 


4.3 


(5.4) 


11.2 (14.6) 


1.6 


(2.0) 


7.2 (9.5) 


2.1 


(2.7) 


9.0 (11.8) 


— 1 

P/P 


0.36 


(0.37) 


0.68 (0.70) 


0.48 


(0.50) 


0.76 (0.78) 


A/A 


0.39 


(0.39) 


0.73 (0.74) 


0.49 


(0.51) 


0.80 (0.81) 


dN B _ B /dy 




16 


.6 


14 


L0 



Table 3: Multiplicities and total transverse energies without and with feed-down contribu- 
tions for 5 % most central Au+Au collision with y / ijv r /v = 130 and 200 GeV. Results with 
feed-down are presented in parentheses. For details, see the text. Also the net-baryon number 
obtained from Eq. @ is shown. 

comes from the weak decays of E, H and Q. (Note that S° decays electromagnetically.) 

As seen in Table the hyperon feed-down contributions remain small for pions and 
kaons 7 but are quite significant for lambdas and especially for protons. The number 
of thermal excitations of heavier hadrons grows with increasing T^ ec and some of the 
decay chains of these heavy hadrons end to proton or lambda increasing their yields 
strongly. Table |3] shows that e.g. the number of A's goes up by a factor of four when 
Tdec is raised from 120 MeV to 150 MeV. A further ~30-40 % increase to p and A 
yields comes from the hyperon feed-down while the ratios p/p and A/ A remain almost 
constant. 

Our pion multiplicities in Table El computed with Td ec = 150 MeV, agree with both 
STAR and PHENIX results in Table El For this Td ec we get, however, a slight excess of 
kaons over the STAR and PHENIX data, as can be seen also in Fig. m For protons, 
our results with feed-down are consistent with the STAR data but the corresponding 
results without feed-down overshoot the PHENIX data. Table 01 shows clearly how 

7 Kaon yields are not affected by the A decays. 
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sensitive the baryon multiplicities are to the decoupling temperature, emphasizing the 
necessity for a high T dec here. 

The baryon multiplicities have also a slight dependence on the number of hadron 
resonance states included in the EoS. It is, however, not clear how many resonance 
states should be taken into account. It depends e.g. on the nature and dynamics of 
the phase transition and on the lifetime of the hadron gas. There is no guarantee 
that the number ratios of hadrons at hadronization are the thermal ones. Since the 
(anti) quark and gluon structure of hadrons varies, formation of some states could be 
favoured over others. In the hadron gas interactions among hadrons would drive the 
densities towards the thermal ones. Since the lifetime of hadron gas is not very long, 
and the temperature drops fast, especially the formation of heavier states might not 
be very effective. 

Even though the baryon multiplicities are very sensitive to T dec , the net-baryon 
number remains independent of T dec , since we have incorporated the conservation of 
net-baryon current in the hydrodynamic equations (JHJ). The final net-baryon number 
thus equals the initial one shown in Table H3 Also the net-proton number is almost 
independent of T dec , the change is only a few per cent when T dec changes from 150 to 120 
MeV. Thus, the measured net-proton number gives a further probe for the computed 
initial state, governed by the saturation scale in our framework. As seen in Tables |2] 
and EJ our results for the net-proton number agree with the STAR and PHENIX data 
very well. 

5.2 From RHIC to LHC: spectra and integrated observables 

The benefit of a closed calculation for the initial state is that we can predict the 
hadron spectra, multiplicities and transverse energies in central A+A collisions also at 
the LHC. The initial conditions are given in Tabled Without a dynamical model for 
decoupling (see e.g. Refs. jllUI II 11] ) , however, we cannot predict the value for the 
decoupling temperature. Based on the RHIC results, we would expect T dec to be in 
the 150 MeV range. To show the sensitivity of our LHC results to T dec , we shall quote 
them with T dec = 120 MeV also. 

In Fig.E3we show the charged particle spectra in 5 % most central Pb+Pb collisions 
at ^snn = 5500 GeV from our hydrodynamical calculation with T dec = 150 . . . 120 MeV 
(the shaded band labeled "Hydro"), as well as from the pQCD+fragmentation with and 
without energy losses. For the pQCD spectra without the energy losses, the shaded 
uncertainty band around the solid line stands for the difference between K\ « 0.65 
from fit I and Ki ~ 1.03 from fit II (see Sec. I4.1j) . In the pQCD calculation with 
the energy losses included, we apply the latter K factor (fit II has a slightly better \ 2 
than fit I), and the large uncertainty band corresponds to the difference between the 
reweighted and non-reweighted results as discussed in the Sec. 14.21 

We see that the hydrodynamical spectra are considerably flatter at the LHC than at 
RHIC. This is due to the stronger flow developed at the LHC. At the same time, the 
suppression in the pQCD spectra is larger than at RHIC 19J. The net result of these 
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Pt [GeV] 

Figure 15: Our prediction for the transverse momentum spectra of charged hadrons at rj = 
in 5 % most central Pb+Pb collisions at the LHC energy y / J/vAr = 5500 GeV. The shaded 
band in the hydrodynamic results shows the sensitivity of our results to T^ ec in the interval 
150. . . 120 MeV. The solid curve labeled "pQCD" corresponds to the pQCD fragmentation 
results where no energy losses have been taken into account. The shaded band related to this 
is the estimated uncertainty in the extrapolation of the K factor to the LHC energy. The 
shaded band labeled "pQCD + E-loss" describes the uncertainty in the pQCD fragmentation 
results with energy losses, see the text for details. 



effects is that the thermal part of the spectrum is expected to dominate over the pQCD 
tail in a px region up to px ~ 4 — 5 GeV that is almost twice as wide as that at RHIC 
where the pQCD tails take over already around px ~ 2 — 3 GeV. This will provide an 
opportunity for more stringent tests of the thermal and hydrodynamical behaviour of 
produced matter, including a wider px range where the elliptic flow should follow the 
hydrodynamical behaviour. 

Figure shows our prediction for the spectra of positive pions, kaons and pro- 
tons at the LHC for the same collision parameters as in the previous figure. The hy- 
drodynamical results are again given with the two different decoupling temperatures, 
T dec = 150 MeV (solid), and T dec = 120 MeV (dotted). We show the pQCD frag- 
mentation results with energy losses using K = 1.03. The LHC results show similar 
systematics between the hydrodynamic and pQCD fragmentation results as at RHIC, 
Figs. El and El the crossover region moves to larger px for more massive particles. 
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Figure 16: Transverse momentum spectra predicted for positive pions, positive kaons and 
protons at y = in 5 % most central Pb+Pb collisions at ^snn = 5500 GeV. The solid 
lines show our hydrodynamic results with T dec = 150 MeV and the dotted lines the results 
with Td C c = 120 MeV. The shaded bands correspond to the pQCD fragmentation results with 
energy losses. 



Integrated observables, total multiplicities and transverse energies, and the multi- 
plicities of different hadron species are given in Table H] for three different decoupling 
temperatures. Even though the slopes of spectra at p? > 1 • ■ ■ 2 GeV depend on T dec , 
the total multiplicity and transverse energy are relatively insensitive on T dec j^J. Like 
at RHIC, only a ~ 6 % change is seen in dN/dy and a few per cent change in dEx/dy 
when changing T doc from 150 to 120 MeV. Note that also at the LHC the pion multi- 
plicity remains quite insensitive to the decoupling temperature. For kaons the change 
in multiplicity is ~ 20 % (see also Fig. HH1 at small pr) but for (anti)protons a factor 
2.8. For pions and kaons the change in multiplicity is negligible if Td ec is raised from 
150 to 160 MeV, whereas the (anti-)proton multiplicity increases ~ 25 %. 

As seen in Tables^andEJ we predict an almost net-baryon free central rapidity region 
at the LHC. The net-baryon number at y = is predicted to be down by a factor 5 
relative to RHIC. Although the proton and antiproton multiplicities are sensitive to 
T dec , also the net-proton multiplicity, dN p /dy — dN p /dy, is practically independent of 
T dec in the range T dec = 160 . . . 120 MeV. Due to the smallness of the net-proton number 
relative to the (anti)proton number, also the antiproton-to-proton ratio p/p = 0.98 for 



33 



rjl [TV T "\ rl 

idee [MeVJ 


1 on 
120 


i rn 

150 


160 


dN/dy(tot) 


4460 


4730 


4840 


dN/drj(XoX) 


4120 


4240 


4290 


dN/dy (char) 
dN / dri{c\iai) 


2760 
2560 


2850 
2570 


2900 
2600 


dE T /dy(tot) [GeV] 


4010 


4070 


4110 


dE T /drj{tot) [CreVJ 


3790 


3710 


3680 


dN* jdy 
dN n ° jdy 


1170 


1120 


1120 


1290 


1240 


1240 


dN n ~ jdy 


1170 


1120 


1120 


dN K+ /dy 


175 


214 


218 


dN K ~ jdy 


175 


214 


218 


dN p /dy 
dN p /dy 


25.8 
24.6 


70.8 
69.6 


88.1 
86.9 


P/P 


0.95 


0.98 


0.99 


dN B _ B /dy 


3.11 



Table 4: Multiplicities and total transverse energies without feed-down contributions 
for 5 % most central Pb+Pb collisions at ^snn = 5500 GeV. Also the net-baryon 
number predicted by Eq. (jHJ) is given. 

Tdec = 150 MeV changes only by a few percent when Td ec = 160 . . . 120 MeV. 

5.3 How independent are hydrodynamic and pQCD spectra? 

Interestingly, as seen in Figs. El and ^1 there is quite a large difference between the 
slopes of the hydrodynamic and pQCD spectra. Our results thus suggest that there 
is a rather narrow px interval in which the dominant hadron production mechanism 
changes from thermal emission at decoupling to independent partonic jets suffering 
energy losses. It should be emphasized that we have not attempted to construct an 
additive two-component model, where the sum of the hydrodynamic and pQCD spectra 
would give the complete spectra, valid over the whole px range. First, in a two- 
component model we should include only the thermalizing partons, the region p sat < 
q T < qrth, in the computation of the hydrodynamic initial state. Correspondingly, 
only the non-thermalizing partons, the region qT > qrth, should be included in the 
pQCD fragmentation production of hadrons. Second, if we consider energy losses of 
jet-producing partons, we should add the energy losses back to the hydrodynamic 
system. Since we do not develop such a two-component model here, and since both 
computations are now initiated with the same partonic cross sections, one should ask 
to what extent our hydrodynamical and pQCD fragmentation spectra in Figs. IDIandfTol 
can be regarded mutually independent, i.e. how much double counting there is if one 
simply adds the two components together. In the regions where one of the mechanisms 
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dominates, this is obviously not an issue but it becomes especially interesting in the 
crossover regions 3 < px < 4 GeV at RHIC and 5 <pr < 6 GeV at the LHC, where both 
components are of equal size. We shall argue below that in these crossover regions 
(and of course at larger values of p?), the hydrodynamic and pQCD spectra are in 
fact biased to different gr-regions of originally produced partons and therefore the 
two contributions now computed are to a good approximation additive without serious 
double counting. 

Consider first the RHIC spectra for 5% most central Au+Au collisions at ^snn = 
200 GeV shown in Fig. El In the pQCD fragmentation calculation without energy 
losses, hadrons at pr ~ 3 GeV are dominantly produced by partons at qx ~ Pt/(z) ~ 
1.1 Pt ~ 5.1 GeV, where (z) m 0.6 follows from the slope of the partonic qx spectrum 
and the shape of the (gluonic) fragmentation functions [T%1 1112j . When energy losses 
are accounted for, the dominating partonic momenta (for a fixed pt) obviously grow 
further, typically by, say 8 , a few GeV. Therefore the hadrons at pt ~ 3 GeV produced 
in pQCD fragmentation + energy loss can be argued to be dominantly produced by 
partons originally at qr ~ 1 GeV. In particular, the small-gr region qx < 3.5 does not 
contribute much. On the other hand, we notice that in the pQCD+saturation com- 
putation of the initial transverse energy, 95% of <j(Et) (see Table comes from the 
region p sat < q T < 3.6 GeV. Were we to exclude the region q T > q Tth = 3.6 GeV from 
the calculation of the initial conditions, the hydrodynamic spectra in Fig.|§]would thus 
decrease only by a few per cent. Similarly, exclusion of the region qT < qrth = 3.6 GeV 
from the pQCD fragmentation + energy loss calculation should not visibly decrease 
the obtained pQCD spectra at pr > 3 GeV, or even at Pt^,^ GeV. Thus we argue 
that to a good approximation the computed hydrodynamic and pQCD fragmentation 
+ energy loss spectra in Fig. |U] at pr > 2 GeV, and especially in the crossover region 
3 1z,Pt ^4 GeV, can be added together without significant double counting. Compari- 
son with the RHIC data also supports this conclusion. 

At the LHC, in the 5 % most central Pb+Pb collisions shown in Fig. 15, hadrons 
at px ~ 5 GeV in the pQCD fragmentation originate dominantly from q^ ~ 8.5 GeV 
partons when energy losses are not accounted for. The average energy losses are now 
larger than at RHIC, due to the denser system formed. Thus, we argue that with the 
energy losses included, the pt = 5 GeV hadrons mainly come from partons which are 
originally produced at, say, qx ~ 11 ... 12 GeV. In particular, the region qx^G GeV 
should not contribute much to the production of 5 GeV hadrons in the pQCD frag- 
mentation + energy loss computation. On the other hand, for the LHC we notice 
that 90 % of the initial ct(Et) comes from p sat < qr < 6.6 GeV and, correspondingly, 
that exclusion of the region qr^Q-Q GeV from the initial state computation would 
reduce the hydrodynamic spectra in Fig. 15 by less than 10 %. Thus, again in the 
crossover region 5 < pt < 6 GeV and at higher px at the LHC, the hydrodynamic and 
pQCD fragmentation + energy loss spectra can be to a good approximation regarded 
mutually independent, and the two components in Fig. 15 added together. 

8 Unfortunately, our current numerical energy loss set up does not allow for a more quantitative 
estimate at pr < qro- 
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The basic assumption in the above discussion is that primary small-gr partons (g*r < 
qrth) thermalize and the large-g^ partons (qx > 9rth) g° through the thermalized 
matter just suffering energy loss but fragmenting independently. We have argued above 
that if we set qxth ~ 3.5 GeV at RHIC and qxth ~ 6 GeV at the LHC, the hydrodynamic 
and pQCD spectra do not significantly change (reduce) from those shown in Figs. EH 
and 15. It should be stressed, however, that the transition from thermalized to jet- 
producing partons does not take place at any well-defined g^th but is gradual. There 
are partons in the region around g^th which are not fully thermalized but are not 
well described by energy-losing independent high-energy partonic jets either. Fig. |U] 
suggests that our hydrodynamic and jet fragmentation components are, when added 
together, not quite enough to fully explain the RHIC data in the crossover region. This 
leaves room for an extra contribution from the non-thermalized partons near qrth (e.g. 
hadronization with the thermal ones) which the present treatment is not sensitive to. 

What makes the final shape of the full hadron spectrum in the crossover region very 
interesting is its dependence on the details of thermalization. A precise experimental 
determination of the shape provides an extra test for the understanding of thermaliza- 
tion of primary partons. 

6 Discussion 

To summarize, we have computed the transverse momentum spectra of hadrons in cen- 
tral and nearly central A+A collisions at RHIC and LHC in a broad pt range. For the 
low-]?T spectra we apply the hydrodynamic framework with pQCD+saturation initial 
conditions [TJJ Ej . From the calculated primary production, supplemented with the 
assumption of adiabatic expansion, we correctly predicted the multiplicities at RHIC 
[171131211. For the high-pr part, we make a separate calculation folding together the 
nuclear PDFs, partonic pQCD cross sections, energy losses and vacuum fragmentation 
functions. 

Assuming thermalization at formation and no initial transverse flow, we have a 
closed framework for computing the initial state for radially expanding boost invariant 
hydrodynamics. Our results show that the main features of the RHIC data on pion, 
kaon, (anti)proton and (anti)lambda spectra, measured in (nearly) central Au+Au 
collisions at p T < 2 . . . 3 GeV, are quite well reproduced with a single high decoupling 
temperature Td ec — 150 MeV. The role of high Td ec for simultaneous chemical and 
kinetic freeze-out is discussed. Also the feed-down contributions from hyperons are 
studied in detail. We emphasize that the net-baryon multiplicity is predicted by our 
calculation and our prediction compares fairly well with the RHIC data. Unlike baryon 
or antibaryon multiplicities, it does not depend on T dec , since the net-baryon number 
is conserved in the hydrodynamical evolution. 

The great benefit of a closed framework is that it can be easily extended to the 
LHC energy. We have presented detailed predictions for hadron spectra in (nearly) 
central Pb+Pb collisions at LHC which depend both on the calculated initial conditions 
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and the hydro dynamic expansion. The overall multiplicity is less dependent than the 
spectra on the details of expansion phase and the predicted net-baryon number as a 
conserved quantity is completely fixed by the calculation of primary production. 

The overall behavior of our hydrodynamic results at RHIC looks somewhat different 
from that at SPS energies. The experimental results for Pb+Pb collisions at SPS are 
reproduced with a decoupling temperature of the order Td ec ~ 120. . .140 MeV [2*2"] . 
the lower values somewhat favoured by the p — p spectrum. That there is a difference 
in the behaviour at SPS and RHIC may not be so surprising since the systems start at 
quite different densities. At RHIC the matter stays in the plasma state longer and as 
the matter enters the hadron phase it flows faster than at SPS. 

One interesting detail, not studied in this work, is how the transverse profile of the 
initial energy density [HO] affects the flow and the determination of T^ ec . As explained 
in Sec. 12 .'2\ here we consider a binary collision profile with a fixed p S at(A \/snn) and 
fixed To = 1/psat- In a central A+A collision, the saturation scales can, however, be 
expected to grow towards the denser region and decrease towards the edges jSH]. This 
would lower the density at center and increase it at the edges. Also, the center of the 
system should form earlier than the regions around it. This may have an effect on 
the transverse profile of initial energy density, as the build-up of transverse collective 
motion can start earlier at center [00] • Furthermore, the form of the initial transverse 
profile should in fact depend also on the cms-energy due to parton production dynamics 
in local saturation [SB]. These studies are left for later work. 

It is quite interesting that with one common decoupling temperature, one can repro- 
duce the gross features of the low-pr RHIC data so well. Certain details, in particular 
the fall of the proton spectrum from hydrodynamic calculation below the data at 
Pt ~ 3 GeV need, however, improvements. As can be estimated in Figs. El and fTT| 
lowering of Td ec to 140 MeV would improve the proton spectrum but it would somewhat 
deteriorate the agreement of the pion and kaon spectra with data. This points to the 
need of separate chemical and kinetic decoupling in describing the details of all spectra 
simultaneously. Studies with separate chemical and kinetic decoupling where the sta- 
ble particle numbers are fixed after chemical freeze-out, show indeed that the spectra 
of pions and kaons become almost independent of the kinetic decoupling temperature 
Tdec,kin, whereas the (anti)proton spectra widen with decreasing Td ec ,kin !"""*J] 172*] . 

In our collinearly factorized LO pQCD calculation of the high-pr hadron spectra with 
absolute normalization, we make use of the K factor's y/s systematics demonstrated 
in [TE]. On top of this baseline calculation, we also add nuclear effects in the PDFs 
[4*7] and the energy losses of hard partons in the framework of eikonalized quenching 
weights [201 02] ■ The amount of energy losses is controlled by the time(path)-averaged 
transport coefficient q. The value of q for central Au+Au collisions at ^/snn = 200 GeV 
was determined in [TO] by fitting the observed relative suppression in the region where 
hadron-type dependence disappears. The current work, supplementing [TO], thus pro- 
vides an important cross-check of the overall normalization of the computed spectra. 

Our pQCD results, extended down to pr ~ 2 . . . 4 GeV at RHIC, show very interest- 
ing relation to the hydrodynamic results. As seen in Figs. "S HTJ1 the pQCD fragmenta- 
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tion + energy loss spectra of charged hadrons and pions merge with the hydrodynamic 
one quite smoothly in the region where the hydrodynamic calculation starts to fall 
below the data. As the slopes in these two spectra are quite different, the region of 
crossing from one mechanism to the other is not sensitive to moderate changes in either 
contribution. As discussed in Sec. 15.31 the two results can be to a good approximation 
considered independent in the crossover region and at higher p?. 

In particular, we can confirm that the difference between the proton data at p? ~ 
3 GeV and the hydrodynamic results for T^ ec = 150 MeV cannot be explained by 
pQCD fragmentation + energy loss results considered here. As discussed above, a 
separate chemical and kinetic freeze-out seems to keep the spectra of kaons and pions 
almost independent of but widens the proton spectrum so that the present data which 
ends at px ~ 3.7 GeV can be described with hydro calculation. It would be extremely 
interesting if the measured px range of the proton spectrum could be extended up to 
6. . . 7 GeV. We would expect the pQCD calculation of jet fragmentation with energy 
loss to reproduce the data in that region and the possible deficit of predicted spectra, 
if any, should show up in the pt interval from 4 to 6 GeV. 

As discussed in detail above, our prediction for the high-p-r spectra at the LHC 
is subject to uncertainties arising from the extrapolation of the K factor and trans- 
port coefficient into a new energy and density regime, and from the normalization of 
the quenching weights. In particular, the uncertainties related to using the quenching 
weights with finite kinematics are quite large up to pt ~ 100 GeV [T9"] . The uncer- 
tainty of the K factor can be reduced by scale-optimized NLO calculations jSH] and by 
applying the latest NLO nPDFs |113j and, ultimately, by comparison with the hadron 
spectra to be measured in p+p and p+A collisions at the LHC. Implementation of 
quenching weights into the NLO pQCD framework will be a challenging task for the 
future. The most urgent need would obviously be the inclusion of finite kinematics in 
the computation of the quenching weights. Also the effects of detailed hydrodynamic 
evolution with transverse expansion of the dense partonic system should be imple- 
mented into the energy loss calculation. 

Acknowledgements 

We thank P. Huovinen, J. Rak, D. Rischke, C. Salgado, K. Tuominen and U. Wiede- 
mann for discussions. This work was financially supported by the Academy of Finland, 
projects no. 50338 and 206024. H. Honkanen gratefully acknowledges the financial sup- 
port from the Wihuri foundation. 

References 

[1] [STAR Collaboration], Phys. Rev. C 70 (2004) 054907 |arXiv:nucl-ex/0407003| ; 
S. S. Adler et al. [PHENIX Colla boration], Phys. Rev. C 71 (2005) 034908 
[Erratum-ibid. C 71 (2005) 049901] |arXiv:nucl-ex/0409015| . 



38 



J. D. Bjorken, Phys. Rev. D 27 (1983) 140. 

K. J. Eskola, P. V. Ruuskanen, S. S. Rasanen and K. Tuominen, Nucl. Phys. A 
696 (2001) 715 |arXiv:hep-ph/0104010| . 

J. Adams et al. [STAR Collaboration], |arXiv:nucl-e~/ 0409033 S. S. Adler 
et al. [PHENIX Co llaboration], Phys. Rev. Lett. 91 (2003) 182301 
|arXiv:nucl-ex/0305013|; B. B. Back et al. [PHOBOS Collaboration], 



'13| 
L2, 



arXiv:nucl- ex/0407012 

P. F. Kolb, |arXiv:nucl-th/0407066| in Proc. Strong and Electroweak Matter 2004, 
Helsinki, June 2004, Eds. K.J Eskola et al, World Scientific 2005. 

J. Adams et al. [STA R Collaboration], Phys. Rev. Lett. 91 (2003) 172302 
|arXiv:nucl-ex/0305015| . 

C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 89 (2002) 202301 
|arXiv:nucl-ex /02060TT| . 

S. S. Adler et al. [PH ENIX Collaboration], Phys. Rev. C 69 (2004) 034910 
|arXiv:nucl-ex /0308006|. 

K. Adcox et al . [PHENIX Collaboration], Phys. Rev. Lett. 88 (2002) 022301 
|arXiv:nucl-ex/0109003| . 

B. B. Back et al. [PH OBOS Collaboration], Phys. Lett. B 578 (2004) 297 
|arXiv:nucl-ex/0302015| . 

I. Arsene et al. [BRAHMS Collaboration], Phys. Rev. Lett. 91 (2003) 072305 
[arXiv:nucl-ex/0307003| . 

I. Vitev, J. Phys. G 30 (2004) S79 1 |arXiv:hep-ph/0403"089| ; K. Filimonov, 
J. Phys. G 30 (200 4) S919 |arXiv:nucl-ex/04Q3060| ; U. A. Wiedemann, 
arXiv:hep-ph/0503119; U. A. Wiedemann, in Proc. of Strong and Electroweak 
Matter 2004, Helsinki June 2004, Eds. K. J. Eskola et al, World Scientific 2005; 
and references therein. 

S. S. Adler et al. [PHE NIX Collaboration], Phys. Rev. Lett. 91 (2003) 072303 
|arXiv:nucl-e x/0306021| ; J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 
91 (2003) 072304 |arXiv:nucl-ex/0306"024| . 

C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 90 (2003) 082302 
|arXiv:nucl-ex/0210033| . 

U. W. Heinz and P. F. Kolb, Nucl. Phys. A 702 (2002) 269 
arXiv:hep-ph/0111075 



arXiv:nucl-ex/0408014 



D. Magestro, J. Phys. G 31 (2005) 265 



39 



[16] H. Von Gersdorff, L. D. McLerran, M. Kataja and P. V. Ruuskanen, Phys. Rev. 
D 34 (1986) 794; ibid. 2755. 

[17] K. J. Eskola, K. Kajantie, P. V. Ruuskanen and K. Tuominen, Nucl. Phys. B 570 
(2000) 379 |arXiv:hep-ph/9909456| . 

[18] K. J. Eskola and H. Honkanen, Nucl. Phys. A 713 (2003) 167 
|arXiv:hep-ph/0205048| . 

[19] K. J. Eskola, H. Honkanen, C. A. Salgado and U. A. Wiedemann, Nucl. Phys. A 
747 (2005) 511 |arXiv:hep-ph/0406319| . 

[20] C. A. Salgado and U . A. Wiedemann, Phys. Rev. D 68 (2003) 014008 
[arXiv:hep-ph/0302184| . 

[21] P. Huovinen, P. V. Ruu skanen and J. Sollfrank, Nucl. Phys. A 650, 227 (1999) 
larXiv:nucl-th/98070"76| . 

[22] J. Sollfrank, P. Huovinen, M. Kataja, P. V. Ruuskanen, M. Prakash and R. Venu- 
gopalan, Phys. Rev. C 55, 392 (1997) larXiv:nucl-th/9607029| . 

[23] J. P. Blaizot and A. H. Mueller, Nucl. Phys. B 289 (1987) 847. 

[24] K. Kajantie, P. V. Landshoff and J. Lindfors, Phys. Rev. Lett. 59 (1987) 2527. 

[25] K. J. Eskola, K. Kajantie and J. Lindfors, Nucl. Phys. B 323 (1989) 37. 

[26] X. N. Wang and M. Gyulassy, Phys. Rev. D 44 (1991) 3501. 

[27] L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rept. 100 (1983) 1. 

[28] V. S. Fadin, E. A. Kuraev and L. N. Lipatov, Phys. Lett. B 60 (1975) 50; I. I. Bal- 
itsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28 (1978) 822 [Yad. Fiz. 28 (1978) 
1597]. 

[29] Y. L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641 [Zh. Eksp. Teor. Fiz. 73 (1977) 
1216]; V. N. Gribov and L. N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438,675 
[Yad. Fiz. 15 (1972) 781, 1218]; G. Altarelli and G. Parisi, Nucl. Phys. B 126 
(1977) 298. 

[30] A. M. Stasto, K. Golec-Biernat and J. Kwiecinski, Phys. Rev. Lett. 86 (2001) 596 
|arXiv:hep-p h/0007192|| ; K. Golec-Bie rnat and M. Wusthoff, Phys. Rev. D 60, 
114023 (1999) |arXiv:hep-ph/9903358l; K. Golec-B iernat and M. Wusthoff, Phys. 
Rev. D 59, 014017 (1999) |arXiv:hep-ph/9807513| . 

[31] A.H. Mueller and J. Qiu, Nucl. Phys. B268 (1986) 427. 

[32] K. J. Eskola, H. Honkanen, V. J. Kolhinen, J. w. Qiu and C. A. Salgado, Nucl. 
Phys. B 660 (2003) 211 |arXiv:hep-ph/0211239| . 



40 



[33] K. J. Eskola, H. Honkanen, V. J. Kolhinen, J. w. Qiu and C. A. Salgado, 
|arXiv:hep-ph/0302185[ Sec. 1.4.2. in CERN Yellow Report CERN-2004-009, 
"Hard Probes in heavy- ion collisions at the LHC", Eds. M. Mangano, H. Satz 
and U. Wiedemann. 

N. Armesto, Eur. Phys. J. C 26 (2002) 35 ||arXiv:hep-ph/0206017| . 
D. Kharzeev and E. Levin, Phys. Lett. B 523 (2001) 79 larXiv:nucl-th70108006 



D. Kharzeev, E. Levin and M. Nardi, arXiv:hep-ph/0111315; D. Kharzeev, 

E. Levin and M. Nardi, Nucl. Phys. A 747 (2005) 609 |arXiv:hep-ph /0408050|. 

K. J. Eskola, K. Kajan tie and P. V. Ruuskanen, Eur. Phys. J. C 1, 627 (1998) 
|arXiv:nucl-th/9705015| . 

K. Morita, S. Muroya, C . Nonaka and T. Hirano, Phys. Rev. C 66 (2002) 054904 
|arXiv:nucl-th/0205040| . 

K. J. Eskola, K. Kajantie, P. V. Ruuskanen and K. Tuominen, Phys. Lett. B 543, 
208 (2002) |arXiv:hep-ph/0204034| . 

L. McLerran and R . Venugopal an, Phys. Rev. D 49 (1994) 2233 
[arXiv:hep-p h/9309289| , ibid. 3352 |arXiv:hep-ph/9311205| . Phys. Rev. D 50 
(1994) 2225 |arXiv:hep-ph/9402335 . 



J. P. Blaizot and F. Gelis, Nucl. Phys. A 750 (2005) 148 |arXiv:hep-ph/0405305 ; 
H. Weigert, in Proc. of Strong and Electroweak Matter 2004, Helsinki June 2004, 
Eds. K. J. E skola et al, World Scien tific 2005; J. Jalilian-Marian, J. Phys. G 30 
(2004) S751 |arXiv:nucl-th/0403077| . 

A. Krasnitz, Y. Nara a nd R. Venugopalan, Phys. Rev. Lett. 87 (2001) 192302 
|arXiv:hep-ph/0108092| . 

K. J. Eskola, Nucl. Phys. A 698, 78 (2002) jarXiv:hep-ph/0104058| . 



T. Lappi, Phys. Rev. C 67, 054903 (2003) [arXiv:hep- ph/0303076 . 

A. Freund, K. Rummukainen, H. Weigert and A. Schafer, Phys. Rev. Lett. 90 
(2003) 222002 |arXiv:hep-ph/0210139| . 

N. Ar mesto, C. A. Salgado an d U. A. Wiedemann, Phys. Rev. Lett. 94, 022002 
(2005) |arXiv:hep-ph/0407018| . 

K. J. Eskola, V. J. Ko lhinen and C. A. Salgado, Eur. Phys. J. C 9 (1999) 61 
|arXiv:hep-ph/9807297| ; K. J. Eskola, V. J. Kolh inen and P. V. Ruuskanen, Nucl. 
Phys. B 535 (1998) 351 larXiv:hep-ph/9802350| . 



41 



[48] K. J. Eskola and K. Tuominen, Phys. Lett. B 489 (2000) 329 
|arXiv:hep-ph/00Q2008|; K. J. Eskol a and K. Tuominen, Phys. Rev. D 63 
(2001) 114006 |arXiv:hep-ph/0010319 . 



[49] K. J. Eskola and K. Kajantie, Z. Phys. C 75 (1997) 515 | arXiv:nucl-th/9 610015 . 

[50] K. J. Eskola, B. Mu ller and X. N. Wang, Phys. Lett. B 374 (1996) 20 
|arXiv:hep-ph/ 9509285|. 



[51] P. V. Ruuskanen, Nucl. Phys. A 702 (2002) 259 |arXiv:nucl-th/ 0111057]; K. J. Es- 



kola, Nucl. Phys. A 702 (2002) 249 |arXiv:hep-ph/0111223| . 

[52] K. J. Eskola, H. Niemi, P. V. Ruuskanen and S. S. Rasanen, Phys. Lett. B 566 
(2003) 187-192 |arXiv:hep-ph/0206230| . 

[53] K. J. Eskola, H. Niemi, P. V. Ruuskanen and S. S. Rasanen, Nucl. Phys. A 715 
(2003) 561 |arXiv:nucl-th/0210005 . 



[54] W. Broniowski and W. Florkowski, Phys. Rev. Lett. 87 (2001) 272302 
|arXiv:nucl-th/01060"50| . 

[55] S. D. Ellis, Z. Kunszt and D. E. Soper, Phys. Rev. Lett. 64 (1990) 2121. 

[56] A. Bohr and B. R. Mottelson, Nuclear structure I, Benjamin, New York, 1969. 

[57] M. Gliick, E. Reya and A. Vogt, Z. Phys. C 67 (1995) 433. 

[58] K. J. Eskola, K. Kaja ntie and K. Tuominen, Phys. Lett. B 497 (2001) 39 
|arXiv:hep-ph/00092"46| . 

P. F. Kolb, U. W. Heinz, P. Huovinen, K. J. Eskola and K. Tuominen, Nucl. Phys. 
A 696 (2001) 197 |arXiv:hep-ph/0103234|. 



[60] P. Huovinen, Nucl. Phys. A 715 (2003) 299 arXiv:nucl-th/02 10024 1. 



[61] K. J. Eskol a, K. Kaja ntie and K. Tuominen, Nucl. Phys. A 700 (2002) 509 
|arXiv:hep-ph/01063"30| . 

[62] P. Arnold, J. Lenag han and G. D. Moore, JHEP 0308 (2003) 002 
|arXiv:hep-ph/0307 325|; P. Arnold, J. Lenaghan, G. D. M oore and L. G. Yaffe, 



Phys. Rev. Lett. 94 (2005) 072302 arXiv:nucl-th/0409068 



[63] A. Rebhan, P. Romatschke and M. Strickland, Phys. Rev. Lett. 94 (2005) 102303 
|arXiv:hep-ph/0412016| . 

[64] R. Baier, A. H. Muelle r, D. Schiff and D. T. Son, Phys. Lett. B 502 (2001) 51 
|arXiv:hep-ph/00092"37| . 



42 



65] H. Niemi et al, in hep-ph/0311131| Sees. 4.6.6.-10. in CERN Yellow Report CERN- 



2004-009, "Hard Probes in heavy-ion collisions at the LHC", Eds. M. Mangano, 
H. Satz and U. Wiedemann. 

[66] I . G. Bearden et al. [BR AHMS Collaboration], Phys. Rev. Lett. 88 (2002) 202301 
|arXiv:nucl-ex/011200l|. 



67] P. F. Kolb and U. Heinz, |arXiv:nucl-th /0305084, in R. C. Hwa (ed.) et al: "Quark 
gluon plasma 3", p. 634-714, World Scientific, Singapore. 



[68] T. Hirano, Phys. Rev. C 65 (2002) 011901 |arXiv:nucl-th /0108004 



[69] T. Hirano and K. Tsuda, Phys. Rev. C 66 (2002) 054905 ||arXiv:nucl-th/0205043 ; 



[70] T. Hirano and Y. Nara, Phys. Rev. C 66 (2002) 041901 |arXiv:hep-ph/0208029 ; 
Phys. Rev. Lett. 91 (200 3) 082301 |arXiv:nucl-th/0301042| ; Phys. Rev~C 
69 (2004) 034908 |a7xTv:nucl-th/0307015| ; Phys. Rev. C 68 (2003) 064902 
|arXiv:nucl-th/0307087| . 

[71] P. Huovinen, P. V. Ruu skanen and S. S. Rasanen, Phys. Lett. B 535 (2002) 109 
|arXiv:nucl-th/0111052| . 

[72] P. F. Kolb and R. Rapp, Phys. Rev. C 67 (2003) 044903 | |arXiv:hep-ph/0210222| . 

[73] K. Rajagopal, Acta Phys. Polon. B 31 (2 000) 3021 [Comments N ucl. Part. Phys. 
A 2 (2002 APCPC,549,95-119.2002) 120] |arXiv:hep-ph/0009058| . 



[74] F. Karsch, J. Phys. G 31 (2005) S633 |arXiv:hep-lat/ 0412038 
[75] F. Karsch and E. Laermann, arXiv:hep-lat/0305025, 

[76] G. M. Welke, R. Venugopalan and M. Prakash, Phys. Lett. B 245 (1990) 137. 
R. Venugopalan and M. Prakash, Nucl. Phys. A 546 (1992) 718. 

[77] R. A. Schneider an d W. Weise, Phys. Rev. C 64 (2001) 055201 
|arXiv:hep-ph/0105242| . 



78] P. Huovinen, |al-Xiv:nucl-th /0505036 



[79] P. V. Ruuskanen, Acta Phys. Polon. B 18 (1987) 551. 

[80] F. Cooper and G. Frye, Phys. Rev. D 10 (1974) 186. 

>1] K. Hagiwara et al. [Particle Data Group Collaboration], Phys. Rev. D 66 (2002) 
010001. 

[82] J. Sollfrank, P. Koch and U. W. Heinz, Z. Phys. C 52 (1991) 593. 



[83] J. Adams et al. [STAR Collaboration] , |arXiv:nucl-ex/0311017] 

43 



84 



85 



86 



89 



90 

91 

92 

93 
94 
95 
96 
97 



98 



99 



E. Schnedermann, J. Sollfrank and U. W. Heinz, Phys. Rev. C 48 (1993) 2462 



arXiv:nucl-th/9307020|. 

F. M. Borzumati and G. Kramer, Z. Phys. C 67 (1995) 137 
|arXiv:hep-ph/9502280| . 

P. Aurenche, M. Fontannaz, J. P. Guillet, B. A. Kniehl and M. Werlen, Eur. Phys. 
J. C 13 (2000) 347 |arXiv:hep-ph/9910252| . 

X. N. Wang, Phys. Rev. C 61 (2000) 064910 |arXiv:nucl-th/9812021 . 



Y. Zhang, G. I. Fai , G. Papp, G. G. Barnafoldi and P. Levai, Phys. Rev. C 65 
(2002) 034903 larXiv:hep-ph/0109233| . 

I Vitev, Phys. Lett B 562 (2003) 36 |arXiv:nucl-th/0302002| ; I. Vitev, 
|arXiv:hep-ph/0212109[ also in Chpt. 2 of CERN Yellow Report CERN-2004-009, 
"Hard Probes in heavy-ion collisions at the LHC", Eds. M. Mangano, H. Satz and 
U. Wiedemann. 

B. A. Kniehl, G. Kramer and B. Potter, Nucl. Phys. B 582 (2000) 514 



arXiv:hep-ph/0010289|. 

H. L. Lai et al [CT EQ Collaboration], Eur. Phys. J. C 12 (2000) 375 
|arXiv:hep-ph /9903282|. 

T. Akesson et al. [Axial Field Spectrometer Collaboration], Nucl. Phys. B 209 
(1982) 309. 

C. Albajar et al. [UAl Collaboration], Nucl. Phys. B 335 (1990) 261. 

G. Bocquet et al., Phys. Lett. B 366 (1996) 434. 

F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 61 (1988) 1819. 

R. Baier, Nucl. Phys. A 715 (2003) 209 | arXiv:hep-ph/0209038| . 

C. A. Salgado and U A. Wiedemann, Phys. Rev. Lett. 89 (2002) 092303 
[arXiv:hep-ph/0204221|; M. Gyula ssy, P. Levai and I. Vitev, Phys. Lett. B 538 
(2002) 282 |arXiv:nucl-th/011207l|; X. N. Wang, Z. Huang and I. Sarcevic, Phys. 
Rev. Lett. 77 (1996) 231 |arXiv:hep-ph/9605213 . 



C. A. Salgado and U. A. Wiedemann, Phys. Rev. Lett. 93 (2004) 042301 
|arXiv:hep-ph/0310079| . 

J. Adams et al [STA R Collaboration], Phys. Rev. Lett. 92 (2004) 112301 
iarXiv:nucl-eI703 10004 1 . 



100] C. Adler et al. [S TAR Collaboration], Phys. Rev. Lett. 89 (2002) 092301 
|arXiv:nucl-ex/02030T6| . 



44 



[101 

[102 

[103 

[104 

[105 
[106 

[107 

[108 
[109 

[110 

[111 
[112 
[113 



K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 88 (2002) 242301 
|arXiv:nucl-ex/0112006| . 



K. Adcox et al. [PHENIX Collaboration], Phys. Rev. Lett. 89 (2002) 092302 
|arXiv:nucl-ex/0204007| . 

S. S. Adler et al. [PH ENIX Collaboration], Phys. Rev. C 69 (2004) 034909 
|arXiv:nucl-ex/0307022| . 



S. S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 91 (2003) 072301 
|arXiv:nucl-ex /0304022 1 . 



B. B. Back et al. [PHOBOS Collaboration] , |arXiv:nuchex70401006 



I. G. Bearden et al. [BRAHMS Collaboration], Phys. Rev. Lett. 94 (2005) 162301 
[arXiv:nucl-ex /0403050 1 . 

I. G. Bearden et al. [BRAHMS Collaboration]. Phys. Rev. Lett. 93 (2001) 102301 
arXiv:nucl-ex/0312023|. 



S. S. Rasanen, Nucl. Phys. A 715, 717 (2003) |arXiv:nucl-th7 0210007 



P. Aurenche, H. Delagrange and P.V. Ruuskanen, in arXiv:hep-ph/0311131, Sec. 
4.7. in CERN Yellow Report CERN-2004-009, "Hard Probes in heavy-ion colli- 
sions at the LHC" , Eds. M. Mangano, H. Satz and U. Wiedemann. 

D. Teaney, J. Lauret and E. V. Shuryak, Nucl. Phys. A 698 (2002) 479 
|arXiv:nucl-th/0104041 . 



D. Teaney, J. Lauret and E. V. Shuryak, arXiv:nucl-th/0110037 
R. Vogt, Phys. Rev. C 70 (2004) 064902. 

D. de Florian and R. Sassot, Phys. Rev. D 69 (2004) 074028 
larXiv:hep-ph/0311227| . 



45 



